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Abstract: In this paper, instead of invoking Dark Energy, we try and fit various cosmo- 
logical observations with a large Gpc scale under-dense region (Void) which is modeled by 
a Lemaitre-Tolnian-Bondi metric that at large distances becomes a homogeneous FLRW 
metric. We improve on previous analyses by allowing for nonzero overall curvature, accu- 
rately computing the distance to the last-scattering surface and the observed scale of the 
Baryon Acoustic peaks, and investigating important effects that could arise from having 
nontrivial Void density profiles. We mainly focus on the WMAP 7-yr data (TT and TE), 
Supernova data (SDSS SN), Hubble constant measurements (HST) and Baryon Acoustic 
Oscillation data (SDSS and LRG). We find that the inclusion of a nonzero overall curva- 
ture drastically improves the goodness of fit of the Void model, bringing it very close to 
that of a homogeneous universe containing Dark Energy, while by varying the profile one 
can increase the value of the local Hubble parameter which has been a challenge for these 
models. We also try to gauge how well our model can fit the large-scale-structure data, 
but a comprehensive analysis will require the knowledge of perturbations on LTB metrics. 
The model is consistent with the CMB dipole if the observer is about 15 Mpc off the centre 
of the Void. Remarkably, such an off-center position may be able to account for the recent 
anomalous measurements of a large bulk flow from kSZ data. Finally we provide several 
analytical approximations in different regimes for the LTB metric, and a numerical module 
for COSMOMC, thus allowing for a MCMC exploration of the full parameter space. 
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1. Introduction 

Our aim in this paper is to check whether spherically symmetric inhomogeneous models 
of our universe can represent a genuine alternative to the Dark Energy paradigm. Dark 
Energy is one of the biggest mysteries of modern cosmology and of fundamental physics in 
general. Despite the fact that it is assumed to constitute about 75% of the energy budget 
of the Universe, we do not know anything yet of its properties except that it dominates 
at late time and that it has an effective equation of state close to —1. We do not know 
if dark energy is only a cosmological constant or if it has dynamical properties, and we 
do not know why it becomes important in the cosmological evolution close to our present 
epoch. A different approach may be taken where we assume that Dark Energy is absent 
or is negligible, and that we do not live at any special time, but rather at a special space 
point. Namely that we live very close to the center of a very large local underdense region 
of the Universe [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. This is of 
course a very radical possibility, which goes against the so-called Copernican Principle, 



which states that we do not occupy any special place in the universe^, and one may argue 
in terms of fine-tunings^ to establish which of the two possibilities is more unlikely; we 
simply take the approach here that observations, rather than postulates about what are 
the structures present in our Universe and what is our location in the Universe, should 
be able to distinguish between the two models. Alternatively, we may think of this as an 
opportunity to test the Copernican Principle on the basis of recent observational data. 

We have in mind a Universe described by a Lemaitre Tolman Bondi (LTB) metric 
(which describes an inhomogeneous region, centered around the observer) matched, at very 
large distances from the observer, to an external Friedmann Lemaitre Robertson Walker 
homogeneous and isotropic metric. We study different LTB profiles, but we always require 
in general the inhomogeneous region to be underdense close to the observer (we call this a 
local Void) . Sometimes there will also be a shell-like overdense structure near the boundary 
of this region (a "compensated" Void). We consider, as a zeroeth order approximation, the 
observer to be at the center, although we do study some aspects of having an off-center 
observer. 

While it is well established that a local Void can mimic an accelerated expansion [3, 
4, 5, 6, 7, 8, 9, 10, 13, 14, 24], whether such a Void can successfully reproduce all current 
cosmological data is still a matter of debate. Most work on the void models have focused on 
reproducing the shape of the ACDM luminosity distance (Dl) versus redshift (z) curve, in 
order to fit the type la Supernova data, but a few studies [14, 15, 16] have included other 
data such as Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations 
scale (BAO). In the present paper we try to perform a much more comprehensive analysis. 
First, we combine fits of several cosmological observations: the CMB measured by the 
WMAP 7-year data, recent Supernova data (the SDSS SN [25], which are a collection of 
SDSS, ESSENCE, SNLS, HST), the Hubble constant measurements (we use here the HST 
values [26]), and the measurements of the BAO scale (BAO) [27]. We also try to include the 
large scale structure data using the SDSS main sample [28] and the Luminous Red Galaxy 
(LRG) subset, DR4 [29], although this required some drastic simplifying assumptions. 

Second, we enlarge the parameter space with respect to previous analyses. So far, 
most fits in the literature have considered a Void embedded exactly in a fiat Einstein de 
Sitter Universe (Eds Void Models) . In [30] voids embedded in a strictly open or flat FLRW 
universe have been considered and fit to the CMB, SN and BAO, although ignoring the 
monopole shift that the CMB experiences in a void (as will be explained in Section 3). Here 
we consider the generic possibility that the outer background FLRW metric has arbitrary 
spatial curvature. Fitting the CMB with an asymptotically curved FLRW model (Curved 
Void Models) allows to fit better the distance to the Last Scattering Surface [31]. For this 
purposes we develop in this paper an analytic treatment of LTB solutions, including the 
derivation of the distance-redshift curve, which allows for curvature in the outer region. 



^Attempts to generalize such a void to a configuration in which the Copernican Principle is not (or less) 
violated, e.g. the so-called Swiss-Cheese universe, so far have failed to agree with observations [19, 20, 13, 
21, 22]. 

^An interesting explanation coining from infiationary cosmology for having such a configuration as a 
generic prediction has been studied in [23]. 



Third, we include the effect of the Void on the redshift of the last-scattering surface, 
which can be encapsulated by an effective change in the CMB temperature monopole. 
This effect is suppressed for small Voids and it was neglected for example in [14], where 
only Voids of size of about 300Mpc//i were considered. However, the effect can become 
significant for larger voids (of order Gpc) and it was taken into consideration in [16]. This 
effect scales as {L/rn)'^-, where L is the size of the Void and rn is the present horizon. 
We also derive the physical quantity that is appropriate for the BAO observations in LTB 
models. This allowed us to accurately compute the physical scales associated with the 
BAO peaks in our inhomogenous model and compare it with observations. 

Finally, we include several variations on the shape of the Void profile. The simplest 
void models essentially have two important physical parameters: the size of the void, 
and the amplitude of density contrast within the void. These models have a relatively 
constant underdensity at the center, followed by a "compensating" overdense shell which 
is then matched to a homogeneous FLRW background. Initially, mainly for the purpose of 
illustration we focus on this simple case. Later however we start exploring more nontrivial 
radial density profiles. Two interesting cases emerge: the non-compensated Voids which 
go only asymptotically to FLRW, and the "compensated" ones which have an additional 
higher underdensity near the center of the void compared to the surrounding region inside 
the void. Both these profiles have a large effect on the monopole of the CMB which can 
accommodate for larger values of the measured Hubble constant and/or locally have a 
higher apparent acceleration providing us with a better fit to the BAO. 

So far, several models have appeared in the literature which mostly differ in the profile 
and the size of the void. It is useful to distinguish between three classes of models: 

I. "Minimal Void" models [14]: These are relatively small in size, extending upto a 
redshift of z ~ 0.1; in this case the main effect is a difference in the expansion 
rate inside (nearby supernovae, z < 0.1) and outside (distant supernovae z > 0.1), 
which can to some extent mimic the effect of acceleration. Since the probability of 
having large voids is exponentially suppressed according to the "standard" analysis 
on growth of structures [32], these models have an advantage over the larger void 
models, but recent studies based on newer supernovae data seem to rule them out. We 
find that allowing for a background spatial curvature makes these models consistent 
with supernovae again, but unfortunately they do not produce a consistent combined 
fit of SN -|- CMB. Thus we will mostly devote our attention to the larger voids. 

II. Large "Compensated" Void models: These typically extend up to very high redshifts 
~ 0.5 — 1.5, and incorporate most of the supernovae; the radial void profile can be 
used to modify the luminosity distance (Di) versus redshift (z) curve all the way up 
to these high redshifts to fit the supernovae data. As we will see, these models can 
indeed be consistent with SN -|- BAO -|- CMB -|- HST. However, in order to have 
a large local Hubble parameter these models require additional features in the Void 
profile. 

HI. Non- Compensated Voids: These match to the background FLRW only asymptoti- 



cally, i.e. it is difficult to say in a sharp way what the Void size is, since they never 
approach exactly FLRW. These profiles do not necessarily have a compensating over- 
dense shell. It turns out that the monopole correction to the CMB temperature in 
such cases helps fit the CMB data significantly. This in turn ensures that much larger 
local values of the Hubble parameter are possible. 

When comparing these models with data the three classes have to be dealt with slightly 
differently. For example, the "Minimal Void" case is well under control analytically using an 
expansion for L/ru <^ 1, which is basically a Newtonian approximation [13, 14]. Moreover, 
all datasets which are located outside this Void (for example the LRG or CMB) can just 
be analyzed within the standard FLRW framework. Finally, when studying effects on the 
Dl — z curve in the outer homogenous region, the corrections can be neglected. 

The Large Void cases (II & III) instead require more care. In particular, if we want 
to employ an analytical approximation, higher order post-newtonian terms in L/rn must 
be included. Besides, the distance curve in the outer (nearly) homogenous region receives 
non-negligible corrections (which have to be taken into account for example when fitting 
the CMB). Finally the data analysis poses non-trivial problems: since many datasets (such 
as the BAO and the large scale-structure) are inside the Void, the usual treatments have 
to be modified accordingly. Finally, for the non-compensated case, since the FLRW metric 
is approached only asymptotically, we have to define what we mean by the background 
FLRW metric. It turns out to be the "effective" FLRW model describing the evolution of 
the last scattering surface. 

While fitting the Void models with the major cosmological data-sets is essential in 
determining the viability of these models, it is also important to find distinguishing features 
of these scenarios which are unique to a spherically symmetric model. In this regard, a 
crucial point is that the observer will in general be displaced from the exact center of the 
LTB metric. This has several observational consequences, such as anisotropic expansion, 
the presence of a large dipole in the CMB, large coherent peculiar velocities ("bulk fiows"), 
etc. We discuss the constraint coming from the CMB dipole on how far away "we" can be 
located from the center of the Void. This is important in estimating the amount of fine- 
tuning involved in such scenarios. The Void also naturally predicts that we will observe an 
"apparent" bulk fiow of all the objects within the void. Furthermore the "bulk velocity" 
should be aligned with the CMB dipole. Remarkably there are now indications of precisely 
such "Dark Flows" of matter [33] around us. It is also worth mentioning that [34] has 
explored a very interesting effect which is present even when the observer is at the centre: 
the kinetic Sunyaev-Zeldovich effect, which can put upper bounds on the size of very large 
Voids. 

When comparing goodness of fit of the different models at stake, we choose to focus 
on the effective quantity Xes = — 21n£, where C is the likelihood. We do not perform a 
Bayesian Evidence (BE) comparison, because the calculation of the BE in the case of the 
Void scenario is not straightforward. We explain this in Section 5.4.1. Secondly, since we 
do find models with a XeS comparable to that of ACDM (at best a difference in x^s of ~ 2 
on 3406 degrees of freedom in the data) but do not find models with exactly the same or 



a lower value of xlsi we do not expect to find significant outcomes of a BE calculation, 
which would anyway still favour ACDM. 

The paper is organized as follows: in section 2 we present the theoretical framework to 
study void models and some general analytical results. In section 3 we discuss qualitatively 
how to compare the void models with data. In section 4 we include a more technical dis- 
cussion, which some of the readers may want to skip, on how to perform such an analysis, 
including the set of parameters that are used in the analysis. In section 5 we describe how 
our numerical code works'^. In section 6 we describe our main results for the simple void 
profiles which help the understanding of how the presence of the void affects the differ- 
ent important observational quantities, and also serve as an illustration of the numerical 
prescription we are following. Next in section 7, we explore more elaborated compensated 
profiles and non-compensated profiles and how it improves the various fits. We also in- 
clude some additional CMB and HST data. In section 8 we comment on other cosmological 
predictions. We conclude in section 9, summarizing our findings, and pointing out some 
unique distinguishing signatures of these models. 

In the appendices we present a number of technical results: in Appendix A analytical 
approximations of LTB metrics and geodesies in different regimes, in Appendix B the 
correction to the CMB temperature in the general case which includes curvature in the 
outer FLRW region and in Appendix C a discussion about the likelihood of multivalued 
functions, which is needed in some part of the parameter space in our Void models. Finally 
in Appendix D we include figures and tables with the results of our analyses. 

2. The LTB model and curvature function 

2.1 The metric 

Traditionally spherically symmetric void models have been studied using LTB metrics 
which are exact solutions of General Relativity that can be studied both analytically and 
numerically. One advantage of using LTB metrics is that the spherical inhomogeneous LTB 
patches can be pasted onto a homogeneous FLRW metric consistently [36, 37, 13]. As far 
as we know, in all analyses that have been performed so far, the LTB metric is patched 
to a flat EdS universe either at a finite "void size" , or asymptotically as we approach the 
Hubble radius. One important modification we are going to implement is to remove such a 
restriction and study solutions where the LTB is embedded inside an open, flat or a closed 
universe. We refer to them as "Curved Void" models as opposed to the traditional "EdS 
Void" models. 

Technically, the LTB metric is described by three different free functions of the radial 
coordinate r, but the picture can be simplified with a gauge choice and a physical require- 
ment. As explained in [37, 13] one function (the "mass function") can be made trivial by a 
redefinition of coordinates, while another function (the "bang" function which sets the time 
of big bang at each value of r) can be set to zero, if we restrict to metrics which become 

^We implemented our code as a module in COSMOMC [35]. We will publicly release our module at 
http://web.physik.rwtli-aaclien.de/download/valkenburg/ 



more and more homogenous in the past. In this case, we are left with only one function 
k{r), the curvature function. Intuitively it is like having an independent scale factor cor- 
responding to each (comoving) radial coordinate, r, which is evolving as an independent 
FLRW metric with a given spatial curvature k{r). The function k{r) is arbitrary and also 
determines the density profile. 

The LTB metric is given in general by 

ds^ = -dt^ + 52 (r, t)dr'^ + R^{r, t){de^ + sin^ Odip^), (2.1) 

Here we have employed comoving coordinates (r, 0, ip) and proper time t. The functions 
S(r,t) and R{r,t) are completely specified by the "curvature" function k{r). S{r,t) can be 
determined from R{r,t): 

Sir,t) = —Ji M , (2.2) 

^l + 2r^k{r)M^ 

where the prime denotes partial differentiation with respect to r, and M is an arbitrary 
mass scale which does not appear in any physically observable quantity. R{r, t) in turn 
is governed by an evolution equation resembling the Hubble equation: as we mentioned 
before, intuitively the LTB model can be understood as a continuous collection of FLRW 
universes, one for each radial coordinate. We may define an r-dependent scale factor and 
a Hubble parameter as: 

R(r,t) , ^^, , d(r,t) 

a(r, t) = ^ ' ^ and H(r, t) = -^^ . 2.3 

Note that this scale factor is the one associated with the angular part of the metric, while 
the radial expansion is in general different. 

We find from the Einstein equations that this scale factor satisfies an r-dependent 
Hubble equation 
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M^ 3k{r)M^ 



a^(r,t) 47ra^(r, t) 
(2.4) 
The above equation is identical to the Hubble equation in a matter dominated FLRW cos- 
mology with spatial curvature, the only difference being the possible r-dependence coming 
from k{r). Accordingly, the solution is also identical to that of a FLRW universe in the 
following implicit form: 

27r 
a(r, t) = —- [ cosh u{r, t) - I] , (2.5) 

t = ^ , [ sinh n(r, t) - u(r, t)] , (2.6) 

3A;(r)3/2M ^ V , ; V , ;j , v ; 

except that the intermediate variable u also depends on the radial coordinate. Note that 
u can be imaginary when k{r) is negative. In this case we have to take the positive root of 
— 1 while taking the square root of k: 

fc3/2 = (Vfc)3 = |fc|3/2^3 ^ _^|^|3/2 fo^ yfe < Q . 



Substituting (2.5) in (2.4) we find that the r-dependent Hubble parameter is given by 

^^_ 3fc(r)3/2M sinhn(r,^) 
^'^" ^/27r [coshn(r,t)-l]2 • ^'■^> 

When the profile is matched exactly to FLRW we have that k{r) = kf, for r > L, where kb 
is a constant and L is the radius of matching with FLRW. So, in the outer region, r > L, 
the r-dependence vanishes in all the quantities and it is possible to see that H(r > L, t) 
corresponds to the "standard" Hubble parameter that is defined for the homogeneous 
FLRW background. In the next section, we will see that the r-dependent Hubble parameter 
at the centre, H(0, t) coincides with the "observational" definition of the "local" Hubble 
parameter using the luminosity-distance vs. redshift relation. We note that our sign 
convention is such that positive k^ corresponds to a spatially open embedding FLRW 
universe, and negative k^ corresponds to a spatially closed embedding FLRW universe. 

2.2 Density Profile and Consistency 

We can compute the matter density and abundance for the LTB metric. Using the exact 
expression for the density function [37, 13], we have 



where Mp is the Planck mass. One can rewrite this as 



p{r,t) = M^ 



3k 1 ^ _ M^ 



27r(coshii — 1) 



a^{r,t) 



(2.9) 



as long k'{r) = 0, which, as we shall soon see, is valid at the centre and in the background 
for r > L. We have defined here the mass scale Mq = M'^M^. From now on we will also 
use the convention that a(-L, ig) = Ij with to referring to the present epoch, so that M is 
fixed by the value of the present energy density p{L,to). 

To have an intuitive understanding of the kind of density profiles that are modeled 
by LTB metrics, let us assume L <^ rn (the Hubble radius) and consider a LTB region 
embedded in an EdS background. Li this case we have 



i^rTTTWTT' where (/,)(t) ^ ^, 
1 -I- [t/to) '^e{r) ottV 



P(^^) « -, , rluJn , ^ ' where {p)(t) ^ -^ , (2.10) 



which is approximately valid as long as 

e(r) = 3k{r) + rk'{r) < 1 and t < to = - ;— • (2.11) 

Mv6vr 

We observe that the EdS behaviour for the density is given by the factor (p) (t) , while 
the fluctuations are provided by the presence of e(r) in the denominator. When e(r) is 
close to its maximum value we have a void, while when it is close to its negative minimum, 
it signals an overdensity. Note that at early times the density contrast 6{r, t) = {p{r, t) — 



(p) (t)) / (p) (t) , defined in the usual way, grows as t"^'^, in agreement witli tlie prediction of 
cosmological perturbation theory. On the other hand at late times, when {t/to)'^''^e{r) ~ 
0(1), the density contrast grows rapidly (and this result is the same as found within the 
Zeldovich approximation ) and diverges. Nevertheless, for our purposes, as long as we do 
not reach this singularity time, LTB metrics adequately capture the effects of non-linear 
stucture formation on photon propagation. 

Now, we are mostly going to be interested in using a k{r) which starts off from a 
maximum at r = and falls off to a constant value at r = L such that 

k'{0) = k'{L) = 0, (2.12) 

HL) = k, = ±-\. (2.13) 

The conditions (2.12) on the derivatives [36, 37, 13] guarantee that there is no density 
cusp at the centre and that there is exact matching to the external FLRW region. Here 
Ofc is the curvature abundance of the FLRW background to which the LTB metric can 
consistently be matched, as explained in section 4.2. (Note that in our conventions fifc > 
corresponds to an open universe.) We will also consider, later on, profiles which are not 
exactly matched to an FLRW. 

The above class of functions describes a spherical void region surrounded by a compen- 
sating shell-like structure. The simplest choice of the curvature function that we employ 
to model the inhomogeneities and fit the supernova data is given by^ 



fc(r) 



i-(£) 



+ kf, for r < L 



(2.14) 
ki) for r > L 



so that the value of "curvature" at the boundary is just kf,, and at the centre 

ko = kb + k^s.^. (2.15) 

We can check that (2.14) satisfies (2.13). There are two important physical parameters, L 
and femax, which correspond to the length-scale and amplitude of fluctuations respectively. 
In the rest of the paper, we will use this profile to gain insight into the basic physics, 
although some of the analytical results are general for any k{r). In section 7, we will 
specialize to more elaborate profiles. 



The exponent of r/L has been chosen to be equal to 4, but the reader may note that any exponent 
n > 1 satisfies the constraints (2.13) and approximately describes the same physical situation of having 
an underdense central core surrounded by a compensating overdense shell. Thus we did not expect n to 
play an important role in fitting the various observations. In general the larger the value of n the flatter 
the void and narrower the structure. For the purpose of illustration we therefore first chose to work with 
(2.14) which was previously shown to provide a good fit to the supernova data [14]. In section 7 we explore 
the effect of changing the profile in greater details. We also note that we could also play with the overall 
exponent in (2.14) which we have set to 2 for similar considerations. 



3. Fitting Observations to Void Models 

3.1 Overview & Basic Physics 

In the previous sections we have discussed the basic properties of the model and in this 
section we want to present an overview of how we proceed to test our model against the 
various cosmological observations. It is well known that the ACDM model provides a 
reasonably good fit to all the important cosmological data. Thus we will use the ACDM 
model as our standard of comparison. We will proceed as follows: we will first look at 
the data coming from SN (Distance vs. redshift), CMB and BAO, together with the HST 
measurement. Next, we will also try to include the Large Scale Structure data, albeit after 
making some simplifying assumptions. These assumptions need to be verified in the future 
and this would require developing a perturbation theory about the LTB metric, which is 
still lacking. 

The main reason why inhomogeneous models can be interesting for the dark energy 
problem is because an underdense region tends to become more and more underdense 
(that is how structure formation works in GR). In other words, the space in an underdense 
region expands a little faster than the rest, and therefore if we happen to be living in a local 
underdensity this effect can reproduce the faster expansion rate as inferred by comparing 
nearby supernovae to high-redshift supernovae. Our first aim is to check, with the advent 
of new supernovae, which inhomogeneous models can still provide a good fit to the SN data. 
In [14] a Small Void profile, extending only up to z ~ 0.1, was proposed but the newer 
supernovae data [25] seem to rule them out. As we will see indeed small EdS void models 
can no longer fit the supernovae data. However a Curved Void models still can, which 
is a first indication that the background curvature may play an important role in Void 
cosmology. It turns out however that the small voids require an open background universe. 
This unfortunately is incompatible with CMB which prefers a fiat/closed universe. Thus, 
we shift our focus mainly to Curved Voids with larger size, and as we will see, they provide 
much better combined fits. 

The second issue we wanted to address is that, although in the literature [14, 15, 16] 
void models can provide a reasonably good combined fit^ of CMB and SN, these fits are 
significantly worse as compared to the ACDM fits. We again find that allowing for a curved 
background improves the fit remarkably by allowing us more freedom on the distance to the 
last scattering surface: compared to the ACDM + curvature model, for the 3403 combined 
data points the EdS Void model has a Ax^ ~ 38, where as for the Curved Void model 
this diff'erence goes down to Ax^ ~ 7. This should not come as a total surprise as it is 
known that the closed universe model, albeit with a very low Hubble parameter, provides 
an excellent fit to the WMAP data [38] . The CMB also receives a non- negligible correction 
to the effective monopole temperature when the void becomes large (this is in contrast 
to the small void scenario studied in [14]), which further changes the distance to the last 
scattering surface. We incorporated this effect in our numerical code to obtain an accurate 
fit to the data. 

''In [14] the p-value (often referred to as goodness of fit) of the combined data sets was found to be 
around 26%. 
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Next, there is the issue of the position of the Baryon Acoustic Oscillations peak in the 
void models. It is well known that the FLRW models without Dark Energy cannot be made 
consistent with the BAO data [39, 14]. However, when the void is large the BAO data fall 
inside the Void, and the distance rulers which are inside the void are very different from the 
usual FLRW ones. We provide a detailed derivation of the relevant distance measurements 
for a general LTB patch. Using these accurate physical quantities, in conjunction with 
having the freedom to choose the background curvature, remarkably, allows us to find a 
good fit to the BAO peak. In fact, for certain inhomogeneous profiles the fit can be made 
even better than the ACDM model, see for instance profile C in table 5. 

In section 6, after discussing the fits to supernovae, CMB and BAO, we look at the 
local Hubble measurement coming from HST. We find that the simplest inhomogeneous 
profiles give us a relatively low Hubble parameter which is difficult to reconcile with HST 
measurements; however we construct other inhomogeneous profiles which provide a higher 
value of h. Depending on which value is taken for HST (different values have been published 
ranging from /i ~ 62 to /i ~ 74), it is possible to provide a consistent combined fit of SN, 
WMAP, BAO and HST. For instance, the combined fit of profile D is only slightly worse 
than the ACDM fit, Ax^ ~ 2 for a total of 3406 data points, when using the HST value of 
62.3 ± 6.3 km s'^ Mpc"^ from Ref. [26]. 

Finally, we discuss compatibility with Large scale structure data coming from the 
SDSS main sample and from the LRG data. The situation here is involved, because in 
void models the matter density is no longer homogeneous and we cannot directly apply 
the known FLRW results on growth of structures. Although progress have been made 
[40, 41, 42] in understanding these issues in LTB models, it is still going to be extremely 
involved to analyze the large scale structure data (LRG and SDSS) in LTB models both 
theoretically as well as numerically. In this paper we took a much easier approximate 
route to deal with this problem: we treat the underdense region in the LTB patch as an 
effective open FLRW model for the purpose of computing the growth of structures relevant 
for SDSS. We discuss our findings and possible caveats in section 6.5. 

In the following subsections we discuss some of the physical quantities that are relevant 
for fitting our model to the observational data, and explain how they can be calculated in 
the framework of LTB metrics. 



3.2 Distance vs. Redshift relation for Supernovae 

The Void models became interesting in the first place because they could fit the SN data 
without Dark Energy. Each supernova is a data-point which measures a luminosity distance 
Dl at a given redshift z. We have included the D^ — z relationship in our MCMC code. 
This firstly involves numerically solving for a radial photon trajectory: 



dt{r) _ R'{r,t{r)) ^ ^^^^^ 



dr 



l + 2/c(r)M2r2 
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Next, the redshift as a function of the radial coordinate is obtained by solving [43, 37] 

"- , '^ + -''^' (3.2) 

dr 



Finally, the angular distances of the supernovae can be calculated quite easily as they are 
simply given by [43, 37] 

DA{r) = R{r, t{r)) and Dl = {I + z^Da • (3.3) 

Thus, once t{r),z{r) is known, so is DL{r). Combining D^^r) and z{r), we obtain the 
usual Dl{z) plot. 

For the supernova data the key physical parameters in our model are kQ,L,ki) and 
to, or equivalently the central density contrast of the void, Sq, the boundary redshift, Zb, 
the background curvature abundance, Q^j and the background Hubble parameter, -ffout- 
(We will discuss in more technical detail the relation between the LTB metric parameters 
and the more physical parameters in section 4.2.) For small voids (zf, < 0.05) and/or 
small inhomogeneities {5o < 0.1), very good analytical approximations for the redshift as a 
function of the radial coordinate exist, see Appendix A.l. We have compared our numerical 
results with the analytical results in Appendix A.l in figure 5. As we can see, both for 
small redshifts and small inhomogeneities, the numerical results agree very well with the 
analytical approximations. 

3.3 Local Measurement of Hubble Parameter 

For local Void models, the essential point is that at the centre of the void the Hubble 
parameter is slightly larger in comparison to the value it has in the external region at the 
same LTB time-slice. In this context we note that the HST measurements go only up to a 
redshift of around 0.1, so that they lie well within the core underdense LTB region. Another 
important remark is that, although we provided an intuitive radially dependent expression 
for the transverse Hubble parameter, observationally the local Hubble parameter, Hq, is 
defined as 

//- = lim^^=5^55_MP^. (3.4) 

In section 4.4 we show that the intuitive definition of the Hubble parameter coincides 
with the formal definition (3.4), i.e. Hq = H{0,t). We also explain how the local Hubble 
parameter is related to //out via 6q and Qk in a general way which is independent of 
the details of the profile. These are therefore the only quantities needed to fit the HST 
measurements. 

The fact that the local Hubble parameter is larger than the background Hubble pa- 
rameter in LTB models is rather useful because, as we will see, the value in the external 
region needed for the CMB fit is very low, and incompatible with HST measurements. It 
turns out that for the simplest LTB profiles the enhancement of the Hubble parameter is 
still not enough to be consistent with HST, but we have found other LTB profiles which 
can yield a much higher local value of the Hubble parameter making it compatible with 
HST, at least using the HST value of 62.3 ± 6.3 km s~^ Mpc~^ from Ref. [26].. This will 
be discussed in section 7. 



12 



3.4 The Cosmic Microwave Background Radiation 

Performing an analysis of the CMB data might seem a challenging taslc in the light of the 
presence of large scale inhomogeneities such as the one we are considering. For the purpose 
of this paper however we are going to assume that we are located very close to the centre, 
so that by symmetry there are no additional anistropic effects on the photon redshift.® So, 
in fact, the CMB can be analyzed in terms of an effective FLRW observer. Let us explain 
this in more details. 

The formation of the CMB spectrum depends primarily on three factors, the primordial 
spectrum of metric fluctuations, the epoch of matter-radiation equality, and the speed of 
the sound waves in the baryon-photon plasma at the epoch of recombination which provide 
us with the characteristic physical length scale Lg , the sound horizon of acoustic oscillations 
at the time of recombination. Throughout this paper we consider an adiabatic primordial 
power spectrum characterized by an amplitude A, a spectral index rig and a running a. 
Both the epoch of matter-radiation equality and the sound velocity at recombination can 
be determined if we know the composition of the universe at the recombination time. Now, 
since the last scattering surface is located in the FLRW background, all we need to know 
are the energy densities of the various components of the universe (i.e. that of curvature, 
dark matter, baryons and radiation) today in the background region. Then we can just 
use the standard FLRW codes to "recreate" the recombination epoch from the background 
universe today. 

There is however one subtlety: While -f/^out, f^fc.out and ilb,out specify the first three of 
the energy density components, for radiation we need to compute the photon temperature in 
the background, Tout, which is not given by the usual value Tq = 2.726i^. The temperature 
seen by the background FLRW observer is in general different, which means that effectively 
the last scattering surface is located at a different redshift. This effect was neglected in [14], 
since the Void considered there was much smaller and the correction to the temperature 
was tiny, since it goes as (L/r/joj,)^- In this paper we calculate the ratio Tout /To which 
boils down to calculating the extra redshift an LTB observer sees compared to the FLRW 
observer. We use numerical results for our MCMC code, but we ensured that our numerics 
agreed with analytical approximations, as discussed in Sections 4 and 5. 

Finally, we need to evolve the physical length scales at the CMB epoch to angular scales 
that we observe in the sky today. This requires the knowledge of the angular distance. Da, 
of the last scattering surface. Fortunately, in the LTB metric this is precisely given by 
R{r,t{r)) which reduces to its FLRW form^ with k{r) = ki, at last scattering. Thus the 



®This symmetry is broken when considering perturbations, and in fact an efTect that we should consider 
is the Integrated Saclis- Wolfe (ISW) effect, due to the secondary effect of density fluctuations at late time. 
Such an effect is absent in EdS, but it is expected to be generically present in a Void model, because 
the growth of perturbations is different, in general. Since we do not have a treatment of the growth of 
perturbation in LTB, we keep the FLRW treatment in the code. The reader should remember, however, 
that the low-^ part of the CMB spectrum cannot be fully trusted because of this approximation. 

^There is actually one non-trivial point about this calculation. Although R{r,t{r)) reduces to its FLRW 
form at last scattering, the function t{r) is different in the LTB patch compared to the FLRW background. 
To account for this mismatch in time, we have to readjust the time coordinate in the effective FLRW model. 
This is explained in more details in section 4. 
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CMB spectrum that we see in the LTB model is identical to the background effective 
FLRW model, and we modified our MCMC code accordingly. To find the parameters for 
this effective FLRW model we need to relate the outer parameters to the inner parameters, 
as a function of the Void profile, density contrast and radius, as we will explain in Section 4. 
Finally, in order to perform fits to the CMB data we also included the other usual pa- 
rameters, which do not depend on the presence of the void. Thus the complete list of param- 
eters in our model relevant for CMB is given by {Ag, Ug, Og, Hout, ^6,out) ^k, ^o, ^b-, t-, ^sz}i 
where r is the optical depth at re-ionization epoch, and Aq,z is the normalization of the 
Sunyaev-Zeldovich template from [44], which is treated as a nuisance parameter since 
from the WMAP observations alone it is not possible to distinguish between primordial 
anisotropics and those induced by the thermal SZ-effect [45]. 

3.5 Large Scale Structure 

The correct way to treat the Large Scale Structure (LSS) data would be to compute the 
growth of perturbations in an LTB model. This is a challenging task, since the per- 
turbations around LTB have not been fully studied, and since they involve a non-trivial 
interplay between radial and transverse modes. Although this has been approached by 
some authors [40, 41, 42] a full treatment is still missing and we do not attempt it here. 
We take the more modest approach of treating perturbations in an effective FLRW model 
to get a sense of how easy or difficult it is going to be to fit the LSS data in the LTB 
models. As mentioned before, in fact, the LTB metric can be qualitatively understood as a 
continuous collection of FLRW universes, one for each r. The transverse Hubble constant 
H{r,t) is given by equation (2.4). Dividing all terms by Hir^to)^, we can define effective 
abundances rim(^) and J^fc(^) via 



SttM^ 



^2(r,to) 



a3(L,to) ^ 3fc(r)a2(L,to) 
a3(r,to) 47ra2(r, to) 



H\r,tQ)a'^{r,to) 
SttM^ 3fco 

^'=(^) = ^^M^WM^' ^'-'^ 

where we are continuing to use the normalization a(L, to) = 1- bi the case of an uncom- 
pensated void that only asymptotically goes to FLRW, we replace a(L,to) with a(r*,to)i 
as described in Appendix B. 

When the void is very large such that zb ^ Zmax, where Zmax is the highest redshift 
for a given LSS survey, and the void profile k{r) is sufficiently fiat near the origin, then 
it seems a reasonable approximation to calculate the matter power spectrum inside the 
void treating it as an effective open FLRW universe. Thus in our MCMC simulations to 
estimate the growth of structures we indeed use the effective open universe model with 
parameters Jlm.in = ^mir = 0), flk,m = ^k{r = 0) and i^obs = ^in = H{0,to)- We 
also assume that Qb/Qcdm remains a constant throughout the universe. Clearly such an 
approach does not take into account the different behavior between radial and transverse 
modes and this should be taken into account in a more refined analysis. 
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To complete the story, apart from the parameters discussed above we also need the 
two other usual "nuisance" parameters to capture the non-linear structure formation: 6^, 
the square of the bias (between luminous and dark matter) in the region outside the void, 
and QnU which parameterizes small scale non-linearities on the matter power spectrum. 

3.6 Baryon Acoustic Peak 

What is measured by several collaborations [46, 27] is a feature in the correlation function 
of the galaxy distribution, corresponding to the sound horizon Ls at CMB time. This 
can be quoted as a measurement of the combination (A^)^A2; at some redshift (so far at 
0.2 and 0.35), where lS.9 is the angle in the sky and Az is the interval in redshift space 
corresponding to the comoving sound horizon Ls- 

Thus, we have to construct what is the value of the quantity (A0)^Az in our LTB 
model. In order to do this, we first have to compute the value of the sound horizon Ls 
corresponding to our cosmological parameters in the relevant location. For this purpose 
we again use the effective abundances (3.6) evaluated at the appropriate r corresponding 
to the location of the BAO observations. Next we compute Ls using a Boltzmann code 
(CAMB) [47] under the approximation of an effective FLRW model. Finally, we calculate 
(A0)^Az using the LTB metric, taking into account properly that the expansion in the 
radial and transverse directions are different. This is explained in detail in the next section. 

3.7 List of the parameters 

To summarize, let us enumerate all the parameters of our model that are needed to fit the 
five observations (SN, CMB, HST, LSS, BAO) we discussed: (i) -ffo.out, this is the Hubble 
parameter outside the void in the FLRW background at time to; which is the time for 
the LTB observer located at r = 0, (ii) J^m.^out, this is the matter (dark matter+baryons) 
abundance in the FLRW background today, (iii) ilfe,out) this is the baryonic abundance in 
the FLRW background today, (iv) (5o, the density contrast at the centre compared to the 
matter density in the background, evaluated today, (v) z^, the redshift corresponding to 
the boundary of the void, (vi) ng, (vii) a^, and (viii) A^, are the spectral tilt, the running, 
and the amplitude of the primordial power-spectrum respectively, (ix) r, this is the optical 
depth at re-ionization epoch. For the matter power spectrum we always marginalize over 
(x) 6^, the square of the bias (between luminous and dark matter) in the region outside the 
void, and over (xi) Qnib^, where Qni parameterizes small scale non-linearities on the matter 
power spectrum. Finally, for the CMB we always marginalize over the Sunyaev-Zeldovich 
template normalization, (xii) A^z- 

4. Technical results: Relating LTB models to observations 

In this section we address some important technical results and approximations which are 
relevant for our analysis. However, due to their slightly technical nature. The reader who 
is interested mainly in the basic picture and results may decide to skip them. 
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4.1 The EflFective FLRW model for CMB 

As explained before, the CMB spectrum obtained in the LTB model is identical to the 
spectrum obtained in an effective FLRW model. 

This is because our observed CMB spectrum primarily depends on four physical in- 
gredients: (A) The primordial spectrum, characterized by amplitude As, spectral index Ug 
and possible running of the spectral index Og ; (B) The epoch of matter-radiation equal- 
ity, which has an important effect on the growth of fluctuations as sub-Hubble density 
fluctuations start to grow as soon as matter starts to dominate. (C) The sound speed at 
the epoch of recombination, which gives us the characteristic length scale for the acoustic 
oscillations, and (D) the angular distance of the last scattering surface which relates the 
angular spectrum that we observe in the sky with the physical length scale of the acoustic 
oscillations at CMB. Now (A) is independent of the "late-time" background cosmology, 
and therefore it makes no difference whether we replace the LTB model with FLRW or 
not. For (B) and (C), since the last scattering surface lies in an FLRW region, all we need 
to know is the energy budget of the FLRW background (the different energy components) . 
We can express this as quantities at a time to = 2/(3iifout) and then back-track and re- 
create the recombination epoch. Finally, in order for the FLRW model to reproduce the 
LTB spectrum the angular distance of the last scattering surface (D) must be identical to 
the LTB model. Since the angular distance in an LTB model is just given by the metric 
function R(r,t{r)), and the function R{r,t) reduces to its FLRW form in the background, 
it seems that we should indeed be able to replace the LTB model with the background 
FLRW model as far as the CMB spectrum is concerned. 

So, we need to find out which FLRW parameters describe an effective observer which 
sees the same CMB sky. Now, R{r, t) has the same functional form as an FLRW model once 
we substitute k{r) = kh in the outer region, so if we want the same angular distances in the 
outer region, we have to find an effective FLRW model such that the trajectory function 
t{r) that appears in the argument of R{r, t) is also identical to the inhomogeneous scenario. 
In general if t(r) coincides in the outer region it will be different in the inner region. So 
the FLRW observer at r = will have a time coordinate t(0) = t^g, different from the LTB 
observer's time to- For similar reasons, the background radiation temperature, which is 
one of the quantities that we need in order to recreate the last scattering surface, is going 
to be different from the observed CMB temperature To, because a photon which passes 
through the inhomogeneous LTB patch suffers an extra redshift compared to its passage 
in an FLRW region. 

The prescription to find t^s and the redshift is as follows. Given an LTB observer 
located at radial coordinate r = and time t = to, we can go backward in time along 
a geodesic t{r) from the centre outwards. Once we cross the Void radius, we reach the 
background FLRW region. We can easily compute the time coordinate, the redshift and 
the angular distance at this boundary, let us denote them by th,Zb and Dj^i, respectively. 
We can now define the effective FLRW model as follows: Let us integrate forward in time in 
the background FLRW model along a geodesic till we reach the center at r = 0. The final 
value of time t(r = 0) = tcff defines our effective fictitious observer. For any object which 
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lies outside the void, this fictitious FLRW observer measures the same angular distance as 
the LTB observer. The time teff will be in general different from to- 

Having defined our effective model and observer, let us return to the question of the 
appropriate background temperature. This is given by the following 

Teff = Tobs:,^^ , (4.1) 

where Tobs = Tq = 2.726 K, and zph is the redshift of the boundary, r = r^, for the fictitious 
FLRW observer. In subsequent subsections we show analytically (for small Zf, and/or small 
(5o) how to compute teffi^eff from Iq^Tq. 

To summarize, the effective FLRW model that we need to use is the same as the 
background FLRW model, but the "present" epoch is no longer given by Iq, but is shifted 
by a small amount to igflf- Accordingly, in our CAMB code we need to rescale the Hubble 
parameter, and the abundances so as to reflect their values at teffj and not at to: 

Ho,cS = H{rb,tes), (4.2) 

^o,out \ / R{rb,to) 



where at r^, we can substitute k{r) = k^,. The radiation temperature is given by 4.1 and we 
keep the ratio il^/ridm constant throughout the universe. 

As it turns out, in practice {^rn,cS-, -f^o,eff} — {^m,out, -ffo,out} for most of the parameter 
space. 

4.2 Relation bet^veen the parameters {-f^o,out, ^m,out, (^Oi 2:;,} and the LTB param- 
eters {fcfo, fcoi ^0,-^} 

In this section we relate the parameters in our void model to the parameters that will be 
used in the Monte-Carlo runs to obtain fits for the different observational data. 

Although the first set of quantities, {i^o.out) ^m,out> (^O) -^b}, are physically more trans- 
parent, technically it is easier to work with parameters characterizing the LTB metric. 
Therefore, first we will explain how the LTB parameters, {kb,kQ,tQ,L}, can be obtained 
from the above "input" parameters. Next, we will show how the different observational 
quantities can be computed from the LTB parameters. 

As we have mentioned, we choose a(L,io) = Ij ^-e. set the value of the scale factor 
today in the background to be one^. Setting a = 1 in (2.9) and (2.4) it is now easy to 
calculate the matter and curvature abundance in the background FLRW cosmology: 



Pm.out = Mg =^ ^-m,out = ^IT , (4.4) 

1 _|_ ■-'"-b 

3fcft 

and Vtk,ont = 1 - ^m,out = — ^^ • (4.5) 

-^ "I" 47r 



This simply corresponds to making a specific choice for the arbitrary mass parameter M: 



3jJ0,out 



B7r(l + 3fc6/47r) 
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Thus knowing O,m,out or fifc^out, we obtain 



kb 



4 Qk 



stt 1 — r^fc 

Further, (2.5) determines Ub in terms of kb'. 

Ub = cosh" I 1 + 



(4.6) 



"ikb 
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From (2.6) we can now obtain ig from the background Hubble parameter: 



(4.7) 



to 



i^o,out V 47r 



(sinh Ub — Ub) 
(coshtif, — 1)3/2 



(4.8) 



It is reassuring to check that as Ub,kb — )• 0, we recover the EdS hmit for the relation between 
the Hubble parameter and proper time. 

We can also obtain ko from the central density contrast. Dividing the central density,(2.9), 
by the background density Mq, we obtain Sq: 



Po — Pout 
Pout 



3kn 



27r(coshiio — 1) 



3kh 



27r(coshtio — 1) 



/^ sinh no — uq 
sinh Ub — Ub 






-1. (4.9) 



where uq = u(0,to) and where we have used (2.6). Since we already know kb,Ub, this 
implicitly determines uq in terms of 6o, which in turn determines /cq via (2.6): 



sinh Uq — uq sinh Ub — Ub 



k: 



3/2 



k 



3/2 



(4.10) 



We have just described an algorithm to obtain {kb, ko, to} from {-ffo,outi ^m,outi <^o} • We can 
also check that using (4.7), (4.8), (4.4), (4.9) and (4.10) we can solve the inverse problem 
as well. 

Finally, we obtain z{L) = zb using a numerical iterative procedure: Since z{r) is only 
obtained after explicit integration of the geodesic equations in a metric for a certain L, we 
employ the following procedure. We have to set some L, and define a new L{L) as the radius 
r for which z{r) = zb- Blueshift may occur in the shell near the edge of the LTB metric, 
such that r{zB) is not unique. But since blueshift will not occur in the embedding FLRW 
metric, it is unambiguous to define L[L) as the largest of all possible r{zB)- Numerical 
inversion is then applied to find the right L for which L[L) = L. In practice this means 
that the code performs about twenty integrations before the right L is found as a function 
of Zb, and L is usually accurate up to one part in 10*^. Since the inversion of L{L) already 
involves solving the geodesic equations, we have the full solution {t{r),z{r)} as soon as 
L{zb) is obtained. This process determines the duration of the calculation, which is of the 
order of 10~^ s and has negligible impact on the duration of the MCMC runs. 
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4.3 The implicit function / and its derivatives 

Although for fitting the various observations with the void models we resort to numerical 
computations, analytical approximations which are valid in different regimes serve as an 
extremely useful check on the accuracy of the numerical codes. In the LTB model (or in 
general in any inhomogeneous model) there are two important parameters, L/r^^or (where 
riior is the horizon radius) determining the size of the inhomogeneity and k{r)^ which is 
related to the spatial curvature/density contrast. In principle, we can perform a perturba- 
tive power-series expansion in either or both the variables. The leading order correction in 
i/^hor is the Newtonian approximation which is a good approximation for small voids. On 
the other hand, usual cosmological perturbation theory in small potentials is closely linked 
to an expansion in k{r) [13]. It turns out that, to obtain these perturbative expansions, 
it is convenient to introduce a function /{k'^'^r'^), where r stands for the EdS "conformal 
time" : 

T = (Mt)5 and 7= [—] . (4.11) 

This function and its derivatives appear on numerous occasions while solving for the pho- 
ton trajectory as well as obtaining relations between different physical quantities. Let us 
therefore try to see why this function is useful in a little more quantitative detail. 
f{x) is implicitly defined as 

l + /(x) = ^(c°'^^^-l) where x=^/2^6(sinh'u-n) . (4.12) 

X 

The function is expandable in a power series of the form: 

°° 1 —3 23 —947 

/(x) = y R^x"", with R2 = — ,Ri= —— , Rq = -— TT-r , R. 



L^ '' ' ^20 2800 ' 504000 ' 388080000 ' 

n=0 

(4.13) 
and so on. What this suggests is that approximately unless x > 20, it will be a very good 
approximation to replace / with just the leading order terms. For future purposes we also 
define its derivatives as 

In terms of this function we can re-express the scale factor as 

a(r,t) = ^7V[l + /(fc(r)TV)] , (4-15) 

To see how useful the function / or rather its Taylor series truncations can be, let us 
compute xo = k^TQ^y'^, which is the maximum value of x that appears in the evolution. 
(Both k{r) and r are largest at the centre.) Let us consider first an EdS background with a 
central underdensity contrast Jo = —0.5. In terms of / we have a rather simple expression 
for the density contrast: 

'^0 = (tT^) - 1 ^^ 1 + /o = (1 + /out)(l + So)-'/^ , (4.16) 
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where 



Now for EdS, font 



2_2\ 



/o = f{koJ Tq) and /out = fihl Tq) . 



/o ^ 0.26 for 60 = -0.5 



(4.17) 



By inverting the function /(x), we find that this corresponds to x ~ 6 and it is clear that 
a truncation of / up to a couple of terms will give rise to an excellent approximation. We 
have verified this numerically. 

Next let us consider the case when we have a curved background, i^k = 0.3 open 
universe. Using (4.5) and (4.7) we then have 



Uh = cosh ( 1 + 



2Qk 



1.23 for Ofe = 0.3 . 



(4.18) 



From (4.12) we can numerically determine that this corresponds to Xf, ~ 1.5 and 1 + /out ~ 
1.08. This gives us /o ~ 0.2 =^ xq ~ 9 which is still sufficiently small. Similarly for 
r^fc = —0.3 closed universe we also find /o ~ 0.2 =^ xq ~ 9. Numerically we find that at 
least up till \i^k\ < 0.3, both for open and closed universes the truncation of / up to 4th 
order terms is a very good approximation. 

4.4 Local Hubble parameter 

The locally observed Hubble parameter is defined via 



Hq ^ = lim -—dL{z) 

^ z-^Q dz 



lim -^dA{z) 
z^o dz 



(4.19) 



where the last equality holds because as z — t- all the different distance measures coincide. 
From the geodesic equations (3.1,3.2), we then have 
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{l + z)R'{r{z),t{z)) 



=0 R'{r{z)Az)) Vl + 2E{r) 



z=0 



lz=0 ' 



(4.20) 



where we used that for z — t- we have k'{r) = by definition, which locally corresponds 
exactly to an FLRW universe, for which indeed R' / R' = R/R. 

For fitting the inhomogeneous LTB models we need to be able to compute the observed 
local Hubble parameter, Hq, in terms of the input parameters of the model. This can be 
done in the following way: First, as discussed in section 4.2, we can obtain the LTB 
parameters k^^ki^ and to from (5o,-f^out and Jl^. Using (4.9) we can then obtain uq. Now, 
using (2.6) and (2.7) we find an expression for Hq in terms of no, to: 



Ha 



1 sinhno(sinhno — no) 



to (cosh no — 1)^ 
This gives us a prescription to calculate the local Hubble parameter in our model. 



(4.21) 
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Alternatively, we can rewrite the local Hubble parameter in terms of /'s: 



rr 2 A + /o + roVfco/i,o A 2 4 + ToVA:o[l + /o] 

""^ - 3to I TTib J - 3^V 4[l + /o]3 ' ('•'') 

where the subscript 0, simply means that / and its derivatives are evaluated at r = 0. This 
expression turns out to be more useful for numerical evaluations. 

4.5 Redshift & The "Background" Temperature 

We have explained before why the background temperature that we have to use in the 
MCMC code to compute the CMB map which is fit to WMAP is different from the observed 
CMB temperature. We also provided a numerical prescription on how to compute the 
effective background temperature, but for a better understanding of the relation between 
To and Tefj , here we provide an analytical treatment of the same in the perturbative regimes. 
Also, since the difference between the observed and the effective temperatures is very small 
for most of the parameter space, the analytical results served as a crucial check on our 
numerics. 

Since we are able to treat the problem analytically in the regime of small k{r) by 
perturbatively expanding the functions involving /{k^'^r'^), we were able to compute the 
correction to the monopole temperature in this regime. For technical simplicity, here we 
will focus on void models embedded in EdS backgrounds, please see appendix B for more 
details. First, let us compare the difference in redshift between an LTB and FLRW observer 
that are both placed at the centre of the void, r = 0, and have the same observational epoch, 

t = to. 

The leading order term in f = Mr ~ r/rn vanishes at the boundary [13], and we have 
to look into O(f^) corrections. The result for an outer EdS metric, derived in Appendix A, 
is given by 



1 + z = {1 + ze) exp 



2ri 
— - + 2aR2j'^TEfk{r) + 2a^R2j^ki{r) 

TE 



(4.23) 



where ze is the redshift for the EdS metric i.e. , with A: = 0. We have also defined the 
following quantities: 

TE = To-af , (4.24) 

^nir) = / dff'^k{r) , and, (4.25) 

n = -a [R2j'^{Tpk + 2aro/ci) - e/ca/S] . (4.26) 

Note that te is simply the conformal time in the EdS metric and that ri is the time 
difference between the LTB metric and the FLRW metric, with the same observer's time 

to- 

The other limit when we can reliably compute the redshift correction is for small f, 
and it is sufficient to keep only the 0{r^) term. This corresponds to Next-to-Newtonian 
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approximation. Again, the results are derived in Appendix A.l and here we quote the 
result: 

l + z={l + ZE)exp{^ f dr rfiri^^k) + —[f{TS'k) +T'^^kh{TS'k)] 

I T6 Jo T 



+ ^ / ^,-,-[3^.^.^;^(^.^.fc) + 2(r^7^fc)V2(T^7^A;)-/(r^7^^^)] 



_2.,2;, ^ ^_2.2,„\ , o/_2.2,,\2^ ^_2.2,„\ f^_2.2, 

LIT r|OT -■'»*■•--'-'' ./_.,;,. -I- .__, .^. *,_-,, 

.2.,2t\ , =r /.,2_2 ^ /_ 2.2; 



X [1 + firh'k) + rk'^^^Mri^^k)]} . (4.27) 

Now, as explained in section 4.1, the effective FLRW observer must have a different 
observer time, t^s, compared to the LTB observer time to to get the angular scales to 
match. In conformal time, the difference is precisely given by 

Tcff = To - n , (4.28) 

since ri was the difference between the LTB and the FLRW patches in conformal time 
that a photon takes to reach r = from r = ri,. This is given by (4.26) in the small k{r) 
approximation, while in the small f approximation we have (Eq. (A. 7) in Appendix A) 

ri = -a (ffir^j^k) + 207^0 f dr rkfiir^j^k) [l + fiT^^k) + fk'^\^h{T^j^k)] 

(4.29) 

Now, in EdS we have 1 + z ~ t~^'^. Therefore, 



1 + ZE = {I + z,s) (^j , (4.30) 

so that finally we can obtain Tefj from Tq using (4.1). 

For small voids, z < 0.2, in EdS background the Next-to-Newtonian approximation is 
reliable, but for larger voids it is no longer valid. On the other hand, for curved FLRW 
backgrounds and density contrasts that we will ultimately be interested in, linear pertur- 
bation theory cannot be trusted. Thus we have to resort to numerical computations in the 
general cases. However, the two approximate limits indeed served as useful checks to the 
numerical results. In Appendix A and figure 5 we show the agreement of the numerical 
results with our approximations. 

4.6 The BAO scale 
4.6.1 FLRW case 

Following [27] the BAO data consists in fitting two numbers, the ratio 

for two values of 2 = 0.2 and z = 0.35. Here Ls is the comoving sound horizon scale at 
recombination and Dy is a combination of angular and radial distance defined as follows: 

Dv = [{l + zfDlD,]'/\ (4.32) 
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In FLRW the "radial distance" is simply given by 



z 



where H{z) is the Hubble rate expressed as a function of z: 

H{z) = Hoy^n^il + z)3 + f^fe(l + z)2 + n^ = Hoh{, 



(4.33) 



(4.34) 



To define the angular diameter distance, D^, let us first define the "comoving distance" as 



r{z) 



dz 



dz 



H{z) Ho J h{z)- 
The angular diameter distance for an open FLRW universe is then given by 

sinh(Fo^/f4r(z)) _ sinh(v^/ 



(4.35) 



Da 



Ho^Ahil + z) 
The expression for Dy then simplifies to 



dz 

Hz) 

HoJn:(i + z) 



(4.36) 



Dv 



Hn 



zsinh^(Vn,/4f) 



1/3 



(4.37) 



In particular a rather useful ratio to consider is 

Dv{0.35) 



n 



Dv{0.2) 



(4.38) 



This quantity only depends on {Qm, ^k, ^a} and therefore provides a rather useful bound. 
We can check that we cannot reproduce the measured ratio without A. In fact even the 
conventional ACDM does not fit the number very well. The measured value [27] is about 
1.812 lb 0.060 , while the ACDM value is 1.67 (with Q\ = 0.75). An open empty universe 
gives about 1.5. 

Having defined how to compute the "BAO observable" in FLRW models, let us now 
describe how to proceed in the more general case of the LTB metric, especially when the 
BAO data is inside the Void. 

4.6.2 Computing the model-independent observable, (A^^Az)^/^ 

Most papers on BAO observations quote numbers for Dy , and this is what we defined above 
for FLRW models. However, what is actually measured by different collaborations ([46], 
[27]), is a model-independent physical observable, the product (A6)'^Az at some redshift 
(0.2 or 0.35). A6 is an angle in the sky, and Az is an interval in redshift corresponding to 
the comoving Ls evaluated at the redshift relevant for the BAO measurements. Only [48] 
gives a value of Az alone and [49] points out that in the future Az (the radial BAO) may 
be more constraining for void scenario, but we do not discuss that here. We are first going 
to present an algorithm to obtain the quantity [A9)'^Az in a general cosmological model. 
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Let us start with the acoustic scale L5. This is a scale for physical processes which 
happen at z ~ 1100. It is usually given as a comoving length and its value is about 
100Mpc//i. 

Given Ls we can reconstruct the value for {/S.9)'^/S.z at some redshift zbao- Let us 
do it first in an FLRW model. This goes as follows: The physical scale at the epoch of 
z = zbao is -^5 ^'^ = 1/(1 + z-Qxo)Lis so that the FLRW angular distance is 

T phys J 

and in terms of the expansion rate ffFLRW ^^ have 



'^S = gPLKW(,,,„) ^ (4.40) 



Therefore 






(4.41) 



which agrees with [50]. 

One can now check that this quantity reduces to 



z 



''L. 



(A^'A-)'/' = ^FLH%- (4-42) 

for FLRW universe. 

So, now if we consider [27], the measured numbers are: Ls/ Dv{^-2) = 0.1905 ± 
0.0061, L5/Dv'(0.35) = 0.1097 ± 0.0036, which is all we need to have the two datapoints 
(A6'^Az)^/3(z = 0.2) and {M'^/\zY/^{z = 0.35). We actually use a 2 x 2 covariance matrix, 
as in [27] in the MCMC code.^ 

Having explained how we can obtain the model independent (A^^Az)^''^ values from 
0's, our next task is to determine a way to compute this quantity for the LTB metric. 

4.6.3 (A6l2Az)i/3 in LTB 

We want to find what our model predicts for {/S.6)'^/S.z at some redshift zbao- So, first of 
all, given a set of cosmological parameters for the outer FLRW region, we can compute, 
using the usual fitting formulae from Ref. [51], the sound horizon in our model Lg '^^^ at 
recombination time, and the redshift to recombination, z^ec- What we really want is the 
sound horizon at the radial coordinate tbao which corresponds to the BAO observations 
(z = 0.35 and z = 0.2). 

We assume here that 2;rec(^) and Lg ''^^'^{r) are not uniform inside the void, but at 
each radius they are approximately determined by the usual fitting formulae from Ref. [51], 



^Another way of expressing this, for instance at z = 0.35, is to consider the measurement _R(0.35) = 
g^ao89) ~ 0.0979 ± 0.0036 [46], where Dm is the comoving angular diameter distance. From WMAP we 
know that the first acoustic peak is detected at an angle of 0.59°, which is 0.01 radians. This means that: 
Ls/ Dm{1Q9>Q) ~ 0.01 , in a model-independent way. Therefore: ^ /f'Jp 3g^ i; '^jsor ^'^'^ ^'^ *^ ^^^ compute 
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except that we have to use Qrn{r), Clk{r) and Clb{r) determined by the equations (3.5, 
3.6). This additional r-dependence is a crucial difference between the LTB and the FLRW 
models, and it is also the main reason why the LTB models can give a better fit to the BAO 
as compared to their FLRW cousins. Moving on, the radius r and time t are determined as 
a function of zbao by the geodesic equation describing the geodesic between the observer 
at (r = 0, to) and the BAO scale at r(2;BAo))^('2^BAo)- 

Then, we just need to see how the scale L^ '^^^ is stretched at the epoch zbao and po- 
sition tbaO) and how this translates to a prediction for AO'^Az. In the transverse direction 
the physical size of the scale at the BAO time is: 

^transverse ^ ^^BAO T^^^rec .^_^^. 



while in the radial direction it is 

^radial ^ ^BAU ^^s'"'"' ' (^-44) 



radial _ -^BAO y-LTB,rec 

R'{r{zBA0),t{Zrec)) ^ 



where we used the short notation R{r{zBAo), t{zBAo)) = Rbao and R'{r{zBAo),t{zBAo)) = 

-^BAO- 

Finally from these two numbers we have to reconstruct the predictions for A^ and Az. 

This is as follows: 

T transverse i 

^^ = ' TTB = ^JTT ¥^^4™'"' ' (4-45) 

Z)LTB i?(r(zBAo),i(^rec)) ^ 

where we have used the fact that: D^^ = Rbao- For the radial direction (in the approx- 
imation E{r) <^ 1) using (3.2) leads to: 

Az = {1 + 2;BAo)-RBAO^^radial , (4.46) 

where Ar^adiai is the coordinate distance which corresponds to the physical length U^'^^^^ 
at (r = rBAO) t = ^BAo)- The two are related by: 

L'^ '"^ {r = rBAO,t = iBAo) 1 ^LTB.roc U A^\ 

Ar-radial = ^, = ^777 TTT ^^g . (4.47 

^BAO R'{r[ZBAO),t{Zrec)) 



Putting everything together we get 

1 



{A9^Azy/^ 



(I-Fzbao)-Rbao^ 



R'irizBAo),t{Zrcc))R'^{r{zBAo),tiZrec)) 



1/3 rLTB 



(1 + Zrcc) 

(4.48) 

In the above equations we need to specify what tree is, i-e. at what coordinate time 
recombination happens at a given position r. We define, for equation (4.48), tree as the time 
for which the volume element in the position r is diluted by a factor (1 -|- Zj.cc)/il + -sbao): 

1 + Zree = (1 + .BAo) ( ^,,^ ^BAogLo \ ' (4.49) 

V R' (tree , "rBAO ) R [tree , rBAO ) J 
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These are the equations we use in our code. More precisely, we define the following 
quantity in the MCMC code: 

Q={M'^/\z/zfl^ , (4.50) 

and we compare with the covariance matrix in [27]. 

Finally, let us make a few remarks about the location of the BAO observations. Let 
us consider a "small" Void {i.e. if it extends only up to z = 0.15 at most): in this case the 
BAO data are all in the FLRW region, therefore we can simply use the FLRW limit of the 
above expressions. 

Using the fact that D^ = ra we get in this case the same expression as (4.41). 

In the more general case in which the Void is larger we take the value of Q{z) in the 
full LTB model at the average redshift [z = 0.2 or z = 0.35). This should be a good 
approximation in the limit of very large voids, that is when the profile does not change too 
rapidly in the range 0.2 < z < 0.4. A more refined treatment could be done in principle for 
more general situations, performing some kind of averaging of this quantity in this redshift 
range, weighted for instance by the number density of observed objects: 

^ SQ{z)n{z)dz 

C^averagc- J„(^)^^ ' i*-^^^ 

where n{z) is a function which describes the number density of observed objects at a given 
redshift. 

5. Numerical Analysis 

5.1 Description of our numerical code 

In order to prevent the mathematical limitations of analytical approximations from bias- 
ing our results, we have written a code that performs a full numerical integration of the 
geodesies in LTB space-times. Our program computes the redshift to any distance. The 
geodesic equations and the analytical expressions for all the necessary derivatives of the 
background functions R{r, t) and S{r, t) have been given in section 2. Analytical approx- 
imations have been used to test the results of our code in all possible limits, as briefly 
discussed in the end of appendix A.l. In particular, we have analytical control in two 
regimes: the approximation for small k[r) (which corresponds to first order in the gravi- 
tational potential, as in the usual perturbation theory around an FLRW metric [13]) and 
the approximation for L <^ rhoD which corresponds to a Newtonian expansion. Our code 
reproduces all analytical results up to high accuracy, and calculates all possible quantities 
to O (lO~^) accuracy. 

As input parameters that specify the cosmological model, we use the set 
{^m,out,f^fc,out,-f^o,out,26,'5o,"'s,as,^s,T, 6^,(5n/,^sz}, as defined in section 3.7. For com- 
parison the ACDM model does not have the z\, and 5q parameters, but it has an additional 
parameter Oai so finally it has the set 

{Qm, ^6, ^A,i -f^o,out, "^s, as, As, T, b^ , Qni, ^sz}i with one parameter less than the Void. We 
will also consider in section 7 more complicated Void profiles with one or two additional 
parameters, which will lead to significantly better fits. 
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In solving the geodesic equations, we can choose different options as the integration 
variable, such as the afhne parameter A, the time t, or the coordinate r. We chose r as our 
integration variable since the metric is analytically defined by k{r), and that there is no 
numerical ambiguity in determining when the photon enters the LTB patch. Note also that 
one should not choose ^; as a time parameter along the geodesic because in a significant 
portion of parameter space the shell of the void is such that the photon experiences a 
blueshift; the over-dense shell contracts faster than the Hubble expansion. This also implies 
an additional complication: for redshifts around the radius of the void, the relation between 
distance and redshift is not unique and typically three different distances correspond to 
the same redshift. The Supernovae observations give us a luminosity distance as a function 
of redshift. Hence, the theory may sometimes predict three different luminosity distances 
corresponding to the same redshift. In Appendix C we construct the likelihood for such a 
situation, and explain why such a likelihood forces us to numerically marginalize over Hq, 
as opposed to doing it analytically, for each fit to the supernovae. Such a marginalization is 
customary, since the Hubble parameter acts as an unknown normalization for the luminosity 
distance - redshift curve from supernovae. 

5.2 Interfacing with CAMB and COSMOMC 

By integrating the geodesic equations we can straightforwardly determine {i^m,m, ^b;nn -f^o.in} 
which are the effective parameters inside the Void, while the parameters {^m,eS, ^b,eS, Ho^cS,Tcs}, 
are used to fit the CMB, as explained in the previous sections. Finally, we numerically 
compute {r2b(rBAo)) ^ml'^BAo)) -f^(^BAo)} necessary for fitting the BAO observations. We 
implemented our numerical code as a module in COSMOMC, which we will publicly re- 
lease at http://web.physik.rwth-aachen.de/download/valkenburg/. For any set of pa- 
rameters, COSMOMC first calls our module, which performs the integration and returns 
{^m,eff) ^fe,eff)-f^o,eff)^CMB,cff}- Then CAMB is run with these effective parameters. If de- 
sired, also {rim,inj ^6,in) -f^o.in} are returned by the module, such that CAMB can be called 
a second time, this time for a calculation of the Large Scale Structure power spectrum 
inside the void. This way we perform a full Monte Carlo Markov Chain (MCMC) analysis 
in order to obtain Bayesian estimates of parameter ranges, and simultaneously the best fit 
X^ that can be achieved for a certain model. 

5.3 Datasets 

For the CMB we use the TT and TE correlation spectra from the WMAP 7-year data 
release [52]. For the supernovae we use the SDSS-II Ist-year SN compilation, where we 
choose the mlcs2k2 lightcurve fitter [25]. For comparison we also include a few runs with 
the SALT-ii lightcurve fitter. The former is supposed to be more cosmology independent and 
more conservative, with larger error bars. We use galaxy power spectrum from the SDSS 
main sample [28] as large scale structure data, although we also include some runs fitting 
to the DR4 LRG power spectrum from SDSS [29]. We use the BAO data as presented in 
Ref. [27]. For the local value of the Hubble constant we use the value 62.3 ± 6.3 km s~^ 
Mpc^^ coming from HST measurements as quoted in Ref. [26] (labeled HST62±6)) although 
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later we also perform runs with higher values of Hq viz., 72 it 8 km s ^ Mpc ^ [53] (labeled 
HST72±8) and 74.2 ± 3.6 km s'^ Mpc"^ [54]^° (labeled HSTo,74±3)- 

5.4 Priors 

We use flat priors on the parameters displayed in Table 1 . Some of the parameters are only 
relevant when fitting to the CMB or LSS, namely {As,ns,as,T}, hence these parameters 
are fixed to a central value when these datasets are not included in the fits. We still allow 
-ffo.out to vary even when fitting to SN only, in order to explicitly test that our numerical 
marginalization over -ffout works as expected. The parameters -ffout, ^b,outh^ and ^dm,outh^ 
are defined in the embedding FLRW universe at time to as explained in section 3. In the 
void scenario, ftk,out = 1 — f^m,out, and in the ACDM scenario f^A = 1 — ^m — ^fe- 

5.4.1 Coordinate divergences 

For certain extreme combinations of 5o and zb, it is very well possible that the metric 
functions S{r,t) and R{r,t) go to zero for coordinates in the massive shell surrounding the 
void. This happens when shells at different radii cross, such that the choice of coordinates 
comoving with the matter is no longer valid. Since the region in parameter space for which 
such situations occur is sufficiently far from the region in which we find good fits, we simply 
reject models for which shell crossing occurs. This implies however, that our prior on zb 
and 6o is not exactly as stated in Table 1, but is bounded by some a priori unknown and 
non-trivial boundary, depending on all cosmological parameters. Therefore we still do have 
a fiat prior on zb and 5o, but we do not know exactly on which volume we take the fiat 
prior. As a consequence it becomes rather complicated to calculate the Bayesian Evidence 
with a code such as multinest [55], which is not compatible with points that are rejected 
by a prior during an assessment of the posterior likelihood. In other words, a Bayesian 
evidence calculation with an a priori unknown prior volume needs further investigation, 
which we leave to future work. 

6. Results with simplest Profiles 

In this section we show results for the simplest Void profile given by Eq. (2.14), which has 
an underdensity and a compensating shell and then is matched exactly to FLRW. In the 
next section we discuss more elaborated density profiles, which will allow to find better 
fits. 

We will fit the data with voids embedded either (1) in an EdS background (EdS Void), 
or (2) in a universe containing dust and curvature (Curved Void). In Table 2 we show, 
for the best fit parameters of each model, a list of XeS = ~21n£, with C denoting the 
likelihood, and differences ^XeS with respect to the best-fit ACDM model. Our ACDM 
model is actually a little different from the "Standard" ACDM model: since we want to 
compare our models with ACDM, to make the comparisons meaningful we wanted the 

^"in this case the quoted value is an effective value, since what is actually fit here is the value of l/dA{z = 
0.04), with dA the angular diameter distance. 
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All scenarios: 



Void: 



ACDM: 



Parameter 



Lower bound Upper bound 



^b,outh'^ 


0.01 


0.03 


^^dmiOut"- 


0.01 


0.2 


i:fo,out [km s-i Mpc-i] 


20 


100 


T 





0.3 


log [I0i%] 


2.5 


3.5 


Us 


0.75 


1.0 


Os 


-0.2 


0.1 



So 



n. 



-0.1 



0.1 



Table 1: Free parameters on which we take a flat prior, and their priors: h = iJo,out/100, the 
optical depth to the last scattering surface r, the amplitude of the primordial spectrum of scalar 
perturbations A^, its tilt Ug and the running of its tilt ag. The parameters i?o,out, ^&,out^^ and 
^dm,out^^ are defined in the embedding FLRW universe at time to as explained in section 3. In the 



void scenario, fi/^^out = 1 — fim,out, and in the ACDM scenario JIa 
pivot scale for the primordial spectrum at /cpivot = 0.05 Mpc/h 



i-n„ 



rih- We chose our 



models to resemble each other as closely as possible. Hence we added two additional 
parameters, the running of the tilt, as, and curvature, fifc, to the Standard ACDM model. 
We chose to perform a full MCMC analysis for each of the displayed combination of 
datasets, in order to show in a most explicit and exact way what the effect of including 
each dataset is on the performance of a given model. Let us go through the table from top 
to bottom, addressing the effects of all datasets. 

6.1 Fitting the Supernovae 

One of the difficulties inhomogeneous models face is that, according to the standard anal- 
ysis of the growth of perturbations starting from the primordial inflationary spectrum, the 
existence of large scale inhomogeneous structures with 0(1) density contrasts is exponen- 
tially suppressed (see for instance [50]). This suppression increases as the size of the void 
increases. In [14] a Minimal Void (MV) model was proposed which could fit the SN and 
CMB data available at that time. The Minimal Void, however, failed to fit the newer 
datasets [25]. In our present study we find that if we include spatial curvature, the best 
fit value for the size of the void is still rather small, with zb = 0.221 with a corresponding 
density contrast of 5o = —0.322 and a global spatial curvature ri^ out = 0.693. The Curved 
Void models give a consistent fit to the SN data, having Ax^ ~ —2 compared to ACDM. 
In comparison, the EdS based void model gives a worse fit than the ACDM, Ax^ ~ +2. 
These values indicate that the SN at higher redshifts can be well fit by a very open FLRW 
universe. The void itself plays a role in fitting the low redshift part of the luminosity 
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Datasets 


#dof 


ACDM + ftk 


EdS Void 


Curved Void 


SN 

SN +HST62±6 

BAO+SN 

BAO+SN +HST62±6 
CMB 

CMB+BAO 

CMB+BAO +HST62±6 
CMB +SN 

CMB +SN +HST62±6 
CMB+BAO+SN 

CMB+BAO+SN +HST62±6 
CMB +SN+SDSS 
CMB +SN+SDSS+HST62±6 


288 

289 

290 

291 

3115 

3117 

3118 

3403 

3404 

3405 

3406 

3423 

3424 


(236.5) 

(236.5) 

(239.0) 

(239.1) 

(3371.2) 

(3372.5) 

(3380.3) 

(3608.9) 

(3608.9) 

(3613.6) 

(3614.1) 

(3627.9) 

(3628.4) 


+2.0 

+2.4 

+3.7 

+5.3 

+35.0 

+39.1 

+32.2 

+37.8 

+38.6 

+40.2 

+40.5 

+49.0 

+50.4 


-2.4 
-2.4 

-2.7 

-2.8 

+4.0 

+7.1 

+3.9 

+6.8 

+12.2 

+6.2 

+10.2 

+19.2 

+28.5 


CMB +LRG+HST62±6 
CMB +SN+LRG+HST62±6 


3425 
3426 


(3638.5) 
(3639.1) 


+30.0 
+30.3 


+7.2 
+12.7 



Table 2: Ax^ for different models against standard ACDM allowing for a nonzero curvature. 
The values in parentheses are the actual — 21n>C = Xoff' i^ ^^^ correct normalization such that 



lim^ 



I £ = 1. For complicated likelihood estimations such as for the CMB, the quantity — 21n£ 



is not x^-distributed, but we still denote it by Xeff for an intuitive interpretation. 

distance diagram. 

On the other hand, the marginalized posterior likelihoods for these parameters are 



ZB 



1.48 



+1.45 
-1.40 



, ^0 



-0.5161°;^^^ and ^fc,out 



0.340 



+0.558 
-1.180 



. The discrepancy between the 



best fit values and the central marginalized likelihood values, tells us that zb is bound 
from above by our prior, not by the data. If the void is very large, say zb ~ 3, which is 
well beyond the highest redshift of the SNe, the SN are practically fit by the effective open 
FLRW universe inside the void. 

We also verified that by including only HST along with SN, the goodness of fit does 
not change at all, indicating that the SN alone do not constrain Hq, and the numerical 
marginalization is implemented correctly. 

6.2 Adding the BAO 

In [14] it was noted that inhomogeneous models based on small Voids cannot be consis- 
tent with the BAO peaks observed at z ~ 0.2 and at z ~ 0.35. We find here that EdS 
based void models, see table 2, have a rather poor fit of the BAO peak position, and this 
effect gets even more pronounced when the IIST62+6 data are included: for a combined 
SN+BAO+HST62+6 best-fit the EdS Void model gives a Ax^ ~ +5. However, this situa- 
tion changes quite dramatically once we allow the background geometry to be curved. In 
fact, the best-fit Curved Void model has Ax^ ~ —3 with respect to ACDM. 

While this is encouraging news, we also have to be cautious before drawing conclusions: 
the improvement in the fit is mainly due to the fact that for a Void in EdS a large Zb is 
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needed to fit the SN, while in our model the freedom in 0,k allows a good fit of the supernova 
even with a lower value of Zb in conjunction with having an open outer FLRW universe (in 
fact we saw above that any value for Zf, is allowed). This allows the MCMC code to find 
a good fit to the BAO with z^ ~ 0.3. Since typically around the overdensity the distance 
scales are very different from the ones in the background or in the Void, the existence of 
an overdense shell at z ~ 0.3 provides enough freedom to find a good fit to BAO. This 
is consistent with our general understanding of the problem of fitting BAO without dark 
energy in FLRW models. Unless, the FLRW universes around z ~ 0.2 and z ~ 0.35 are 
significantly different from each other we cannot hope to get a good fit of the ratio of 
the Dy's (4.38) at these two redshifts. Fortunately, LTB models has enough latitude to 
reconcile the two measurements. 

However, it is clear that the good fit that we obtained is a somewhat fine-tuned result: 
first, because it is quite dependent on the shape of the profile of the overdense shell at 
the BAO location, second because if the BAO scale changes rapidly it is no longer a good 
approximation to take the distance scale at one single redshift point, which is what we do. 
More importantly, as we will discuss below, when we add the CMB data, the best-fit model 
will no longer correspond to Zf, ~ 0.3, essentially because the open outer FLRW model does 
not provide a good fit to the CMB. Thus we will not reproduce this very good fit of the 
BAO. As a general lesson however, we see that we can fit the BAO scale well, and that this 
fit becomes especially good if we allow the LTB metric to have some kind of feature at the 
redshift of the BAO. We also saw that Q^ can be a crucial parameter in void models, and 
this will become even more apparent in the next subsection. 

6.3 Adding the CMB 

Apart from the ACDM model, also an FLRW closed universe (without dark energy) does 
fit the CMB data well, albeit with a very small Hubble parameter. LTB models offer a 
unique opportunity to vary curvature of the universe in such a way that the last scattering 
surface be situated in a closed FLRW background, while the local universe mimics an open 
universe. This also increases the value of the local expansion rate and helps us fit the 
supernova data in a manner similar to the EdS based inhomogeneous models. Note also 
that a Minimal Void does not lead to a good fit, because an open curvature in the outer 
region would be needed to fit the Supernova which lie outside the Void, while a closed 
curvature is required to fit the CMB. We find that Q^ plays a decisive role while fitting 
the combined data sets of SN, CMB and BAO. As we can see from Table 2, the O^ = 
best fit model has a Ax^ ~ 40 as compared with ACDM, while our Curved Void model 
gives Ax^ ~ 6. The combined data set contains 3405 data points, and thus if we consider 
the x^ psr degree of freedom, the difference between ACDM and the Curved Void model 
is not large. 

We should also point out that while fitting the CMB, the Void plays an important 
role by affecting the monopole of the CMB temperature: in fact, as explained in section 
3.6, we describe the CMB via an effective FLRW model with a correction on the monopole 
temperature. This effect was neglected in [14], because the Void considered there was much 
smaller and the correction to the monopole was tiny, since it goes as (L/rhor)'^, but when 
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we consider larger Voids this effect is not negligible anymore (see [16]). This effect is such 
that the photons in the effective FLRW have a slightly higher temperature than what is 
observed in the Void, as shown in Table 11. This means that effectively the surface of last 
scattering is closer than in a ACDM universe with the same parameters otherwise. A same 
size of the sound horizon at decoupling would hence be observed at a larger angle for an 
observer in the Void. Note that this correction to the monopole temperature goes in the 
direction of making the CMB fit better {i.e. makes the distance to last scattering larger) 
but note also that this happens only if the Void boundary is roughly 2: < 2. For larger 
Voids, in fact, the correction changes sign, and therefore too large Voids ar disfavoured 
in the MCMC analysis, as can be seen from Figure 7. For a discussion of the sign of the 
correction see subsection 7.1. 

To summarize, the freedom in O^, 60, z-q and ifo,out allows the large Curved Void model 
to fit the CMB almost as well as ACDM, while the EdS Voids cannot be made consistent 
with SN, CMB, and BAO all at once. This is one of the main positive results of our study. 
The overall trend, when following the Table further down, is that ACDM does a better job 
everywhere when the CMB is included, but the Curved Void follows closely. 

6.4 Adding HST, Parameter constraints from CMB+BAO+HST+SN 

Let us now turn our attention to the HST measurement, and in particular see whether our 
model can be consistent with the combined data sets of CMB+BAO+HST+SN. In this 
subsection we discuss the HST value 62.3 ± 6.3 km s~^ Mpc~^ from Ref. [26]. The higher 
values of 72 ± 8 km s"^ Mpc"^ from [53] and 74.2 ± 3.6 km s"^ Mpc"^ from [54] will be 
discussed later in section 7 and will be clearly more difficult to fit. This is actually the 
largest choice of datasets where we can trust the various assumptions and approximations 
we made in our analysis. For instance, we have correctly computed the distance to the last 
scattering surface and the energy content of the universe at last scattering, which makes 
our fit to the CMB reliable. The distance measure to the SN is calculated as good as 
exactly. The value of the Hubble parameter is measured at such low redshifts in HST, that 
the model prediction at r = is a good estimate. And finally the BAO is a relatively model 
independent measurement of a combination of the expansion of space in angular and radial 
(redshift) directions. The only assumption that went into our calculation and prediction of 
the BAO inside the void, is that the local universe at the drag epoch is correctly modeled 
by an FLRW-metric. Thus we can be quite confident about our analysis of the combined 
data set consisting of CMB, BAO, SN and HST62±6- We will discuss in the next subsection 
why LSS data, possibly the only other important observational set, cannot accurately be 
analyzed within our theoretical framework of the LTB universe, although we will try to 
perform an estimate. 

In Figure 7 we show the marginalized one-dimensional posterior probabilities of all free 
parameters on which we had a flat prior except for ^k,out^ which is only a free parameter 
in ACDM and a derived parameter in the Void scenarios, and r2dm,out^^) which is a derived 
parameter for the EdS- Void model. A comparison between ACDM, the Eds- Void model 
and the Curved Void model is provided. In Figure 8 we show the same for all derived 
parameters. In Table 3 we give a breakdown of the total \^ in terms of the x^ against the 
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Model 


CMB BAG SN 


total x^ 


Model 


CMB BAG SN 


HSTg2±6 


total x^ 


ACDM 


3371.5 2.3 239.7 


3613.5 


ACDM 


3371.1 3.1 239.5 


0.4 


3614.1 


Curved Void 


3377.1 3.6 238.9 


3619.6 


Curved Void 


3377.4 4.0 238.9 


4.1 


3624.4 


EdS Void 


3408.8 4.3 240.7 


3653.8 


EdS Void 


3408.7 4.2 240.9 


0.8 


3654.6 


Model 


CMB SN SDSS 


total x^ 


Model 


CMB SN SDSS 


HST62±6 


total x^ 


ACDM 


3371.2 237.1 19.6 


3627.9 


ACDM 


3371.3 236.7 19.6 


0.8 


3628.4 


Curved Void 


3377.0 255.8 14.3 


3647.1 


Curved Void 


3379.3 257.8 12.7 


7.1 


3656.9 


EdS Void 


3411.4 252.8 12.7 


3676.9 


EdS Void 


3411.2 252.4 12.9 


2.3 


3678.8 



Table 3: A breakdown of the various contributions to the total best fit y^ for the different models, 
when simultaneously fit to CMB+BAO+SN (top left), CMB+BAO+SN+HSTgaie (top right), 
CMB+SN+SDSS (bottom left) and CMB+SN+SDSS+HST62±6 (bottom right). When we compare 
from left to right, we see that the inclusion of HST62±6 hardly affects the goodness of fit to SN, 
BAO or SDSS, but it does weigh in on the goodness of fit to the CMB in the Curved Void case. 
Comparing from top to bottom, we see that exchanging the BAO for the SDSS has repercussions 
on the goodness of fit to the SN as well. That is, the (oversimplified) assumption of taking the 
cosmological parameters at ?■ = to fit an FLRW cosmology to the observed SDSS, has a strong 
influence on the allowed range of (5o, handicapping the luminosity distance curve. 



different datasets. In Table 11 we show the best fit parameters and marginalized parameter 
ranges for all three models. 

It is clear that our Curved Void model now provides a significantly worse fit compared 
to the ACDM model: the effect of including HST62±6 is large, even though IIST62±6 only 
includes one datapoint. In comparison the effect of adding the BAO and SN to the CMB 
on the Ax^ of the Curved Void model versus ACDM is not so large; only ~ +3.1 and 
~ +2.8 respectively for 2 and 288 extra data points respectively, please see Table 2. 

Allowing Vtk to vary, significantly opens up parameter space in the sense that this 
model can reconcile the two different constraints on the distance to the last-scattering 
surface and the constraint on Qm^eff^^- However, since the best fit point happens to be a 
slightly closed universe, the much better fit to the CMB is achieved at the cost of a lower h 
and therefore a worse fit to HST62±6) as we see in Table 3. This is exactly the same problem 
that the closed FLRW model encounters in reconciling CMB with other observations. The 
profile of the curvature inside the void as well as the density along the past light cone are 
shown in Figure 1. The inclusion of O^ clearly allows for much larger density contrasts, 
favoured by the SN and CMB. Even after boosting the Hubble expansion rate with the 
inhomogeneous void, the Hubble parameter remains inconsistent with HST62±6- However, 
as we will see in the next section, once we allow the void profiles to have some additional 
features, or deal with very large inhomogeneous models, the problem of a low Hubble value 
ameliorates considerably. Figure 7 shows that the redshift of the boundary of the void is not 
really constrained from above, except from the effect on the CMB monopole temperature 
as discussed in section 3.4, reflected by a thick tail to the right in the posterior distribution 

Oi ZB- 

In Figure 8, one might be intrigued to see that in EdS, the void actually seems to do a 
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better job at fitting HST62±6 (-ffo.m) than the Curved Void. What is going on is a struggle 
between two well-constrained characteristics of the CMB: the position of the peaks and the 
matter to radiation density ratio. The former is most relevant for constraining Q^nfish'^, 
the latter for constraining a different combination of -ffo.eff and Om.eff- Since in this model 
^k,eff = and J^m,eflf = 1, ^m,eff^^ and Hq are strictly related to each other. The matter 
to radiation density ratio favours a value of about rim.eff^^ ~ 0.14, while the position of 
the peaks favours h ~ 0.45. These two requirements are in contradiction with each other 
and eventually, the larger value of -ffo,out and -ffo,in is only a fortunate consequence of an 
otherwise bad fit to the CMB. 

6.5 Large Scale Structure 

For the large scale structure data (LSS) we use both the SDSS main sample, which have a 
mean redshift of z ~ 0.1 and the LRG data which go much further away, with an average 
around z ~ 0.4. For the SDSS data in all the relevant cases, the size of the void, zb, is 
so large that all the datapoints are contained inside the void. For the LRG this is only 
partially true, but we proceed in the same way just to get an estimate. This means that 
we can try to fit the LSS with the approximate effective FLRW universe that is built up 
from the parameters fij^in and ii/^o,m at the centre of the void, as explained in sect. 3.5. 
Let us emphasize the caveats to this approach. To begin with, the large scale structure 
data does not lie at the centre of the void but instead at a slightly higher mean redshift, 
hence our prescription is not so accurate, especially for the LRG. More importantly, the 
growth of perturbations inside an LTB metric is not yet understood. This means that the 
X^ values when fitting the void to SDSS are to be taken only as a rough indication. In 
Table 2 we show the x^ for the different SDSS and LRG runs. The void does badly when 
fitting against the LSS, be it with Ofc.out = or with ^k,ont ¥" 0- In fact, both matter 
power spectra favour a value for the combination ^dm,mh of about 0.2, while the Curved 
Void from the fit with CMB+BAO+SN+HST62±6 has a best fit value of about 0.09. Note, 
however, that the fit is significantly better for LRG than for SDSS. In fact, the void does 
predict too much power on large scales and the LRG data prefer this, compared to the 
SDSS data. 

In Figures 9, 10 and 12 and Table 12 we show the parameter constraints resulting 
from a fit to CMB+SN+SDSS+HST62±6- The value of Qdm,inh is now twice as much, 
with respect to the fit without the SDSS. This is due to the fact that the LSS data prefer 
a higher value of ^dm,inh, which means a lower density contrast. This poses difficulties 
for the SN as well as for the CMB. However, as long as no perturbation theory in the 
LTB-metric has been developed, there is not much to say about the goodness of this fit. 
The strongest conclusion we can draw from this, is that it gives a hint that the void may 
not be in agreement with observed LSS perturbations, but a correct treatment of the 
perturbations is really needed before we can say anything definitevely. The main message 
of this paragraph hence is that such perturbation theory will be of crucial importance for 
testing the Void versus ACDM. 

The effect of the large scale structure on the redshift of the boundary of the void, 
zb, is that smaller radii no longer give a better fit than larger radii, hence pushing zb up 
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towards the regions in which it is not constrained by the data but only by the prior (and 
not even by the monopole effect, since this is smaller due to a smaller density contrast). 
This is apparent if we compare Figure 9 with Figure 7. 

6.6 Baryon density and BBN 

We note that the baryon density in the outer FLRW ri^^eff^^ is constrained by the CMB to 
be in the right range for the Big Bang Nucleosynthesis constraints, and from Table 11 we 
indeed find this number to be very close to the ACDM result. This number tells us what 
the baryon-to-photon number density ratio, 7] = Ufj/n^, was outside the Void. Therefore 
we are constraining t] in the outer region (rjout) and our value of Q,b,csh'^ ~ 0.02 is in good 
agreement with BBN constraint. In this context we note that measurements of the abun- 
dances of the light elements relevant for BBN are taken at high redshift {z ^ 1) [56, 57, 58], 
and therefore r/out is the appropriate quantity to consider. The only measurements which 
are closer to us are the ones relevant to ^Li, which are known to be in disagreement [59] 
with the BBN, so we shall not consider them. Actually [60] we could even take the point 
of view that those measurements are in disagreement with the rest because they are taken 
inside the local Void and we could try and explain the ^Li problem as a consequence of 
the hypothesis that we are living in a Void. 

6.7 Significance of AXgg 

So far we have only quoted the effective ^Xcs = ~2 InC It is still difficult to translate these 
numbers into a fair model comparison. We could not do a Bayesian Evidence calculation, 
as explained earlier. In order to get a taste of the meaning of the quoted x^ values, let us 
quote (approximate) p- values, that is the probability that the goodness of fit can be worse 
than the one obtained, given that the hypothesis is true. This has to be taken with caution 
here: we already stressed that we are quoting an effective x^ > because the likelihood of the 
CMB is not a true x"^ distribution. Nevertheless, we calculate the p-value as if everything 
were exactly x^i for simplicity. The numbers serve to give the reader a feeling for the 
difference between the models. For true p-values one would have to simulate all data and 
explicitly probe how often a worse fit is obtained, when these models are the true models. 
The p-values against the observations are quoted in Table 4. We take only the best 
fit parameters obtained from the fit to CMB + BAO + SN + HSTe2±6; and calculate the 
p-values of these models when fitting to some datasets alone and some combinations of 
datasets, as an illustration. As done in the WMAP paper [61], we consider the CMB as 
three independent tests of the model, namely the high / TT C; spectrum, the high / TE 
Ci spectrum, and the combined low / pixel based likelihood analysis (CMB all low /). As 
could be expected from the x^-values in Table 2, ACDM has the best performing p-values. 
It is still interesting to see that when we demand a rejection at 99% C.L., the Curved 
Void model along with ACDM is not rejected, as opposed to the EdS Void model which is 
rejected by the larger combination of datasets. 
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Model 


CMB TT high / 


SN 


CMB TT high I 


CMB TT high I 


(Bestfit to CMB+ 






+ SN 


+ CMB aU low I 


BAO+SN+HST62±6) 








+ SN + BAG +HST62±6 


ACDM 


5.4 X 10-2 


0.96 


0.29 


4.2 X 10-2 


Curved Void 


5.1 X 10-2 


0.96 


0.28 


2.5 X 10-2 


EdS Void 


2.2 X 10-2 


0.95 


0.17 


1.3 X 10-2 



Table 4: The roughly approximated p-values of the x^ of the best fit models obtained when fitting 
to CMB+BAO+SN+HST62±6: against the different datasets, considering x^-values of the datasets 
alone and in combination. As done in the WMAP paper [61], we consider the CMB as three 
independent tests of the model, namely the high I TT C; spectrum, the high / TE C; spectrum, 
and the combined low I pixel based likelihood analysis (CMB all low /). The p- value we quote here 
for the high I TT spectrum is a factor three lower than the value that the WMAP team obtained, 
probably due to our flawed assumption that the likelihood is actually x^- The last column considers 
all datasets except for the CMB TE spectrum. If we demand a p- value of at least 5%, all models 
are rejected at 95% C.L. by the observations (last column). The high p-values against the SN, 
indicate that the error bars in the MLCS2k2 pipeline from the SDSS SN survey are most likely 
overestimated. These values serve as an illustration of the meaning of the x^ values in Table 2, but 
are by no means to be considered correct. 

7. Results with Profiles with higher Hq 

It is important to realize that with infinite freedom in the function k{r), we can construct 
almost any expansion history on the past light cone^^. Therefore it may be possible to 
fix the problem that we have found in the previous sections, namely to find profiles which 
can fit HST as well as the other datasets. There are two ways this can be achieved. The 
first way is to modify the profile at very low redshift, say z < 0.03, since the observed 
value of h is measured at these low redshifts, and such a change should not affect any 
other measurements, such as SN and BAO. The other way is to try and get a much larger 
correction to the monopole of the CMB. Since it is the CMB which constrains if out to 
be so small, by altering the distance to the last scattering surface significantly, we can 
hope to obtain a different value of if out- This can be done if we use non-compensated 
profiles which asymptotes very slowly to FLRW, and only at high redshifts, say z ^ 3. 
This ensures that the compensating overdense shell, which in general cancels the monopole 
shift from the underdense void, is absent and we in turn can get a much larger monopole 
correction. Further, this modification does not affect other measurements such as SN and 
BAO. We explore in detail the two options in the following subsections. It may be possible 
to combine them together, but we leave that as a future excercise.^^ 

7.1 Modified profiles 

In this section we will consider curvature profiles that exhibit a larger curvature very close 



^^With the exception that one always gets deceleration locally at the origin r = [62, 63]. 

^^Eventually we could perform a principal component analysis, describing the void profile by nodes. This 
way we would find which features in a void profile are most important and best constrained, as is similarly 
done for the equation of state of Dark Energy [64] . 
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Figure 1: Curvature profiles as a function of coordinate radius r (top) and density profile as a 
function of redshift z along a geodesic (bottom), comparing the best-fit CMB+BAO+SN+HST62±6 
to CMB+SN+SDSS+HST62±6- The density is normaUzed to the FLRW density that would be 
measured at that redshift by an observer in a pure FLRW at the same time of observation as the 
observer in the centre of the void. Comparing the figures from left to right, it is clear what the 
restrictive power of the large scale structure inside the void is. 



to the centre of the void, in other words with a more peaked density contrast. Effectively 
this scenario looks like a 'Void inside a Void'. 



We have tested the profiles given below in Equations (7.1-7.4), where profile A is the 
profile considered in the previous sections, given by Eq. (2.14). Profiles B, C, D and E add 
two more degrees of freedom. For B, C and D these are the curvature k^a.x,2 in the inner 
void and its radius L2 = f L. Instead of L2 there is a parameter a in profile E, which 
determines the rate at which the profile converges to FLRW at high redshift. 



The profiles C, D and E only reach FLRW asymptotically. In the case of profile C 
and D this happens at such a high rate, that there is no significant effect, compared to the 
exactly compensated profiles. Profile E on the other hand, only reaches the asymptotic 
metric at a very slow slow rate, and it does so with a very large 'overdense' region, which 
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is discussed in the next subsection. 



Profile A: k{r) = < 
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The profiles for these configurations are illustrated in Figure 2, using their best fit 
parameters. The best fit parameters and the marginalized likelihoods are presented in 
Table 13. We find that the 'Void inside a Void' does indeed improve the goodness of fit to 
the HST observation, but not only that. It also improves the goodness of fit to the BAO 
by a large factor. In fact, the improvement for the BAO is such that it actually makes 
the Void truly competitive with ACDM, as we see in Table 5 and Figure 3. ACDM has a 
best fit x^ value of 3614.1 when fit to CMB+BAO+SN+HST, where profile D is capable 
of getting a x^ of 3616.1. This difference is due to the fact that the profile changes rapidly 
around the redshift of the BAO scale. However, as we said in Section 6, we stress that we 
have to be cautious in this case: when the BAO scale changes rapidly it is not anymore a 
good approximation to take the distance scale at one single redshift point. 

Another interesting feature in these profiles is apparent in Figure 3. The inner void 
clearly has its radius, and hence its shell, somewhere in the middle of the observed su- 
pernovae. Voids are preferred to be large since the inclusion of the SDSS Supernovae [25] 
which filled the gap at intermediate redshifts. But the lesson to learn here, is that the pref- 
erence for large Voids is not due to the shell no longer being allowed inside the observer 
SNe, but only due to the overall shape of the distance-z relationship, fi{z), which favours 
a very open universe up to high redshifts. A feature in fi{z) is not ruled out at all. 
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Figure 2: The best-fit curvature profiles (top), density profiles (second from top) and redshift 
histories (third from top) as a function of coordinate radius r/L, and the density as a function of 
redshift (bottom). The horizontal black lines in the third panel indicate the redshift radius zb of 
the void. In the second panel it seems that for profile E the asymptotic curvature is not reached 
before the surface of last scattering, when we look along coordinate radius r. In the bottom panel, 
however, as a function of the physically relevant redshift, it is clear that also for this profile the 
asymptote is as good as reached at the surface of last scattering. 



We also analyze a class of profiles ^(r) that does not match exactly to FLRW, but 
which goes only asymptotically to FLRW at very large r, i.e. /c(r) — )• /c^ at r — )• oo. It has 
been shown [16] that within this class of profiles we can have a large effect on the monopole 
temperature, due to the fact that such profiles do not have a perfectly compensating shell, 
which in the usual case almost exactly cancels the effect of the Void. Instead, such models 
have an overcompensating shell, such that the total mass contained in the LTB patch is 
actually larger than what would be contained in an equally large background FLRW patch. 

However, if the profile goes to a constant very quickly after some distance L, we get 
an effect very similar to the previously analyzed case. Therefore, to get something non- 
trivial we need a profile which goes slowly to a constant, as exemplified by profile E. It is 
possible to show, already from the analytical approximations given in Appendix A, that 
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Model 




CMB 


BAO 


SN 


HST62±6 


total x^ 


ACDM 




3371.1 


3.1 


239.5 


0.4 


3614.1 


Profile A 


(Curved Void) 


3377.4 


4.0 


238.9 


4.1 


3624.4 


Profile B 




3377.0 


0.2 


237.9 


2.2 


3617.3 


Profile C 




3376.9 


0.7 


237.7 


1.9 


3617.2 


Profile D 




3377.5 


3.6 


233.7 


1.3 


3616.1 


Profile E 




3380.2 


3.3 


241.4 


0.8 


3625.7 



Table 5: A breakdown of the total x^ for each dataset, for fitting simultaneously to CMB + BAO 
+ SN + HST62±6- The main improvement achieved with the 'Void inside a Void' is on the x^ of 
the BAO, more than on that of HST62±6- 
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Figure 3: The best-fit theoretical prediction of a Void with profile C compared to ACDM, for the 
BAO (left) and SN (right). Clearly, this Void profile has less difficulty fitting the BAO than ACDM 
has. An interesting feature in the right figure, is that the inner void has its radius in the middle of 
the SN, which is no problem for getting a very good x^- This implies that the reason that single 
Voids need to be large, is not the disagreement with a feature in the middle of the redshift range, 
but the overall shape of the SN curve that favours a very open universe. 



the dominant contribution to the monopole effect (in the flat case /cf, 
is given by the following integral: 



ko 



k{r)r dr . 



0) for such a Void 



(7.6) 



It turns out that this integral enters with a negative sign in the final expression of the 
monopole, so if k2 is positive this makes the effective CMB temperature lower, rather 
than higher. This is precisely what happens in the compensated case as well; since this 
term starts dominating when the Void has a very large radius, the fit becomes worse (see 
sect. 6. 5). In the case of an asymptotic profile the way to get a positive correction to the 
monopole in (7.6) is to have a k{r) which goes negative and then asymptotes to zero from 
below. In the curved case, the profile has to become more negative than k^, and then 
approach k^ asymptotically from below. 

We find that in such cases the correction to the monopole is much larger than the 
previous 0(1%), going up to 30% — 40%. This drastically changes the other parameters, 
most notably i^o- In figure 4, we see that profile E changes the temperature of the CMB so 
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Figure 4: The 2D marginalized posterior likelihoods of the locally observed Hubble constant 
and the temperature of the CMB for the effective observer, for all profiles: A (solid black), B 
(dashed blue), C (dotted magenta), D (dashed green) and E (dashed red). Only profile E predicts 
significantly different values, notably a higher observed value for Hq. 

much, that the effective outer universe allows for a higher value of //out (and accordingly 
also Hq) compared to the other profiles. 

7.2 Other data sets 

7.2.1 A higher value for H^ 

For the sake of completeness, we performed a brief analysis including more CMB data, and 
other values for Hq. In the previous sections, we used the value Rq = 62.3 ± 6.3 km s~^ 
Mpc~^ from Ref. [26]. The results that we obtained using Hq = 72± 8 km s~^ Mpc~^ and 
Hq = 74.2 lb 3.6 km s~^ Mpc~^, as quoted in Ref. [53] and Ref. [54] respectively, are shown 
in table 6 and table 7. 

Profile E manages to fit the higher HST value Hq = 72 it 8 km s~^ Mpc~^ at 1.7a, at 
the cost of a worse fit to the BAO and SN than profile C does. The change in goodness 
of fit to Hq, is mainly due to the shape of the profile at very large distances. Due to the 
large effect on the monopole of the CMB, profile E allows for a higher overall value of Hq, 
hence also a higher value inside the void. Since the distances at which this monopole-effect 
comes to play are different from the distances at which the profiles fit the BAO and SN, 
it may be possible to combine the good features of profile C and E into one new profile, 
thereby fitting all datasets better at the same time. 

When fitting Hq = 74.2 it 3.6 km s~^ Mpc~^ from Ref. [54], one does not actually fit 
Hq, but rather the value of z/dAiz) at z = 0.04, with dA the angular diameter distance. 
Fitting this quantity is more model independent, since Hq = huiz^Q z/dA, and the actual 
observation is at nonzero redshifts. So the value of /i = 74 makes sense only when assuming 
ACDM in deriving Hq. That this is in fact the case, is refiected in the resulting x^- For 
Profile E for example, the value of i/o,in for both fits (Tables 6 and 7) is -ffo,in = 58 km s~^ 
Mpc~^. But fitting against the efi^ective Hq = 74.2 it 3.6 km s~^ Mpc~^, its x^ is slightly 
less then (74.2 - 58)V(3.6)2 = 20.3, namely x^ = 18.5. 

Putting this minor effect aside, the higher value for Hq with smaller error bars leads 
to no surprises: the Void Profiles assessed here perform in a similar manner overall, but 
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are punished for the lack of a high value for -ffo,m- 



Model 



CMB BAO 



SN HST72±8 total x^ 



ACDM 


3371.8 


2.3 


239.5 


0.3 


3613.9 


Profile C 


3377.1 


4.0 


235.9 


3.6 


3620.6 


Profile E 


3382.0 


5.9 


242.2 


3.0 


3633.1 



Table 6: A breakdown of the total x^ for each dataset, for fitting simultaneously to CMB + BAO 
+ SN + HST72±8, where HST72±8 indicates that we used Fq = 72±8 km s'^ Mpc"! from Ref. [53], 
keeping all other data sets equal to the previous sections. Profile E fits i^o at VS-Ocr = 1.7(t, which 
is good compared to the other profiles, due to its high-density asymptotic profile as discussed above. 
Profile C gives a much better fit to the BAO and SN, such that its overall fit is still slightly better 
than profile E. 



Model 



CMB BAO 



SN HST74±4 total x^ 



ACDM 


3372.7 


1.8 


239.7 


2.1 


3616.3 


Profile C 


3380.7 


1.9 


238.4 


20.7 


3641.7 


Profile E 


3380.2 


3.2 


242.0 


18.5 


3643.9 



Table 7: A breakdown of the total y^ for each dataset, for fitting simultaneously to CMB + BAO 
+ SN + HST74±4, where HST74±4 indicates that we used effectively 74.2 ±3.6 km s^^ Mpc~^ from 
Ref. [54], keeping all other data sets equal to the previous sections. The quantity that is acutally 
fit in stead of Hq is z/dA{z) with z = 0.04, which is more model independent. 



Model 



ACDM 
Profile C 
Profile E 



CMBe BAO 



SN HST72±8 



3531.4 
3537.7 
3537.9 



2.3 

4.2 
2.4 



239.7 
236.1 

242.0 



0.3 
3.3 
6.7 



total x^ 



3773.7 
3781.3 
3789.0 



Table 8: A breakdown of the total y^ for each dataset, for fitting simultaneously to CMBe + BAO 
+ SN + HST72±8j where HST72±8 again indicates that we used Hq = 72 ± 8 km s~^ Mpc~^, and 
CMBe indicates that we use the extra data sets as discussed in the text. 

Again, let us quote the approximate p- values we obtain by assuming that the goodness 
of fit is actually ^^-distributed, and that our Xeff values are the actual x^'s- The values are 
listed in Table 9. In this case, when using 74.2 it 3.6 km s~^ Mpc~^ from Ref. [54], we see 
that the extra freedom in Profiles C and E allows these profiles to obtain similar p-values 
as profile A in Table 4, in spite of the much worse fit to the higher i?o-value. In the last 
column we provide a comparison with the HST62±6-i'uns for these models, where the x? 
is very close to that of ACDM. The higher number of free parameters in the void profiles, 
compared to ACDM, still pushes the p-value down relatively to ACDM. 

7.2.2 Small scale CMB data and additional polarization data 

On top of the higher value for Hq, we also added CMB data from BICEP [65], ACBAR [66] 
and CBIPol [67] along with the already used WMAP data. The resulting best fit x? values 



42 



Model 


CMB TT high I 


SN 


CMB TT high / 


CMB TT high / 


CMB TT high I 


(Bestfit to CMB+ 






+ SN 


+ CMB all low / 


+ CMB all low I 


BAO+SN+ 








+ SN + BAG 


+ SN + BAG 


HST74±4) 








+HST74±4 


+HST62±6 


ACDM 


5.2 X 10-2 


0.96 


0.28 


4.2 X 10-2 


4.2 X 10-2 


Profile C 


4.2 X 10-2 


0.95 


0.26 


1.3 X 10-2 


3.4 X 10-2 


Profile E 


4.1 X 10-2 


0.93 


0.23 


9.7 X 10-3 


2.2 X 10-2 



Table 9: The roughly approximated p- values of the x^ of the best fit models, against the different 
datasets, as in Table 4. The CMB, BAO and SN give no surprises with respect to the previous 
section. In the second column from the right we see that ACDM has no difficulty fitting the higher 
Ho-value, and that the extra freedom in profiles C and E is necessary to keep them at the same 
level with respect to ACDM, when compared to profile A in Table 4. In the last column wc show 
for comparison the p- values obtained in the runs with IISTg2±6, using the best fit parameters from 
those runs. There the difference with ACDM is much smaller, although the void profiles are still 
punished for their extra parameters, since the number of degrees of freedom used for the p-valuc 
calculation is A'dof = A^data points — ^froc parameters- Notc that these p-valucs are very rough estimates, 
which serve as an illustration only. 



are quoted in Table 8. Comparing to Table 6, we immediately see that the inclusion of 
other CMB data has no role in the difference in x^ between ACDM and the Void-models. 
An improvement of the goodness of fit must be searched in the prediction for the observed 
value of Hq, as well as on the large scales of the CMB, which is exactly where the late-time 
ISW-effect is important. However, given the ignorance about perturbation theory in the 
LTB metric, at this point we cannot say precisely what is the impact of the ISW effect, 
and how the void profile should change in order to improve the goodness of fit to the CMB. 



8. CMB Dipole and the "Dark Flow" 

When considering a non-homogeneous large scale model for our Universe, it is unlikely that 
we should occupy a special position, the centre. However the more off-centre we place the 
observer, the more we spoil isotropy. Therefore the observer's position is constrained by 
observations. The most constraining observation is given by a velocity that a comoving 
observer acquires with respect to the CMB frame. Barring special cancellations (the ob- 
server may have a small-scale peculiar velocity which cancels the large scale LTB velocity) 
this puts a strong constraint [9]. So, let us consider our observer to be located slightly 
off-centre, at r = ro- In this case the non-zero radial velocity of the observer will lead to 
a dipole moment in the CMB [9, 14] of the order of: 



6T 



vo = do - dp , 



i.l) 
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where the proper radial distance, do, of the observer is given by 

rro ff 

do = dr 

Jo 



yrT2(iQr)2%) 



and dp is the proper radial distance of an observer at the same position but at rest in the 
FLRW background, such that do = dp but do ^ dp. 

Now, in our profile k{r) remains almost a constant for almost the entire underdense 
region. Assuming we are living in this "constant" underdense region, we have 

, . r° . 1 , ,tanh~H^~^\/2A^^o) ,. , . 

do = aoit) / dr ^ = = ao{t) ; jL ^^ « ao{t)ro . (8.2) 

•JO ,h M9.(M'r\'ih(r-\ My/2ko 



1 + 2(Mr)2/c(r) 

where ao(^) refers to the scale factor at r = 0, and the simplification occurs because u and 
hence R' becomes only a function of time, and Mro is expected to be very small. Taking 
the time derivative we simply find 

do = do{to)ro = Hodo ■ (8.3) 

In the FLRW background, we similarly have 

dp = dpitoyp = Hout do (8.4) 

Thus at the present epoch the estimated dipole moment comes out to be 

^T 27r(coshno - l){Ho - H^niVo , . 

— ^ [Hq- Houtjdo ~ ^ • (8.5) 

We note that uq can be obtained by inverting (2.7). In Table 10 we compare approx- 
ation (8.5) to numerical results for the maximum distance of the observer to the centre 
of the void, for which the dipole is at most equal to the dipole observed in the CMB, 
3.355 mK [68]. For the profiles A and B, which behave similarly at the centre and have 
close to constant k{r) for small r, approximation (8.5) is good up to a few percent. The 
observer is thus constrained to be within 0(10 Mpc) of the centre of the void. 

While the CMB dipole moment constrains our model, the fact that we are located 
slightly off the centre of the void can also account for a rather surprising observation, 
the "dark flow". In [33] the authors observed a systematic difference in the kSZ effect 
from light sources as one varied the angle keeping the distances fixed. The authors had 
already subtracted the contribution coming from the motion of the observer (this gives 
the CMB dipole moment discussed above) and therefore this additional Doppler shift was 
interpreted by the authors to indicate a "dark flow" of all the light sources towards a 
particular direction. This can be understood as follows: If all the objects have a common 
peculiar velocity, say vs, along a particular direction, then the expected Doppler shift will 
be given by 

Azksz ^vscosO , (8.6) 

^'^For an exact calculation of the dipole moment see [?] for instance. For the purpose of estimation 
however it is sufficient to calculate the dipole moment from radial velocities. 
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Profile 


ZB 


L [Mpc] 


?^3.355mK [Mpc] 


r/L 


Eq. (8 


.5) [mK] 


A 


1.024 


4540.964 


17.104 


0.00377 




3.02155 


B 


2.776 


7400.249 


13.814 


0.00187 




3.13626 


C 


2.912 


7586.902 


12.663 


0.00167 




3.22672 


D 


1.012 


4536.971 


5.476 


0.00121 




3.27521 


E 


2.129 


5077.986 


13.394 


0.00264 




3.29353 



Table 10: Dipole constraint on the position of the observer for the best fit parameters of the 
different curvature profiles. rs.assmK indicates the numerically obtained radius at which the dipole 
due to the off-centre position is (up to numerical precision) exactly equal to the observed dipole 
in the CMB, 3.355 mK [68]. The last column gives the dipole as estimated in Eq. (8.5), which 
assumes a constant kir) for small r. 

where 9 is the angle between the position and the velocity vector of a given light source. 
In fact, in [33] a cosine variation of the kSZ effect was extracted on a scale of at least 300 
Mpc//i, and the authors estimated the bulk velocity to be 600 — 1000 km/s from their 
observations, which is way too large to be consistent with ACDM. 

In the LTB model however, the above effect would naturally arise: if we are located 
slightly off the centre of the void, then the light sources which are approximately the same 
distance away from us, will not be at the same distances from the centre. Now, in the void 
region all the objects move away from the centre towards the spherical structure at the edge 
of the void, and more importantly, their outward velocity is proportional to their distances 
from the centre. Thus we, the off-centre observers, will start seeing a cosine modulation 
in their velocities (Doppler shift). Let us try to provide a simple estimate of the effect. 
Approximately using simple Eucleadian geometry we find 

rs^r'sU + '^cosd] , (8.7) 

where rs,r'g are the radial distances of a given light source from the centre and the off- 
centre observer, is the angle of the sources as measured from the observer. Now, if we 
assume that both the observer and the light sources are located in the core of the void 
region where k{r) ss /cq, then we can use the same Doppler shift formula (8.5) that we have 
derived above: 

■ 2'k{Hq- Hont){cos\iUQ-l)rs 2tt{Ho - Hont)icoshuo - l)r'g f tq „ 

''^ = '^s'^ ^, ^ ^, ( 1 + ;rcos^ 

(8.8) 
Please note that since tq ^ rs, the radial directions subtended by the light sources to the 
centre and the observer are approximately the same. 

The term multiplying one is simply a monopole term which is filtered out of the 
kSZ dipole observations, but the second term, containing cosO, provides an additional 
contribution to the kSZ dipole. Thus in the void model we will have 

AzkSZ~cos6' — = cos 61 — . (8.9) 

•^'^0 V ^ / cmb-dip 



45 



Thus we actually see a connection between the dark flow velocity and the dipole measured 
in CMB in an LTB model, which gives precisely a flow of order 10~^ times the speed of 
light on very large scales, roughly consistent with the measurement of [33]. In fact [33] 
mentions that the direction of the bulk velocity also seems to be roughly aligned with the 
dipole velocity, again what we would expect in our model. A prediction of our model is 
that we should continue to observe the same "dark flow" velocity all the way out to Gpc 
scales till we reach the end of the underdense region. 

9. Conclusions 

In the present paper we have investigated the possibility of fitting most of the available 
cosmological data without Dark Energy, with a non-homogenous density distribution in 
which we are located at the centre of a Gpc scale Void that can mimic cosmic acceleration. 
The new features of this analysis are a combined MCMC fit of many datasets (WMAP, 
Supernova, BAO data, HST and an estimate of the Large Scale Structure data) for which 
we wrote a publicly released module in COSMOMC. We also enlarged the parameter space of 
Void models by allowing for the presence of a background curvature, and we have included 
the additional redshift that the photons experience when traveling in a Gpc scale Void 
which modifies the distance to the last-scattering surface. We have shown that the SN, 
BAO and CMB can be fit well already with the simplest Void profile, which includes a 
Void and an outer overdense shell with two parameters, matched exactly to FLRW. In this 
case the best-fit Void extends upto z ~ 1, it has a density contrast of 5 ~ —0.65 and the 
external background FLRW is slightly closed, with fi^ ~ —0.2. The primordial spectra in 
the best fit are consistent with the usual almost scale-invariant ones. Finally, the Baryon 
abundance computed in the outer region is consistent with BBN. 

However, for the simplest profile, the measured value of the local Hubble constant of 
about h ~ 0.6 ~ 0.7 is difficult to fit, since the CMB forces the Hubble constant to be 
very low {h ~ 0.45), in order to fit the angular distance to the last-scattering surface. 
We have then shown that modifying the Void profile by adding more parameters we could 
accommodate for a larger value of Hq, upto h ~ 0.57, in the process obtaining an overall 
better fit to all the cosmological data that becomes comparable to the ACDM fit. This 
can be achieved in two ways. Either one can modify the profile at very small redshift, 
thereby increasing the local density contrast and obtaining a higher local Hubble constant. 
Incidentally, this leads also to a better fit of the BAO data. Or, one can modify the profile at 
very high z introducing a large overdense region which goes to FLRW only asymptotically 
at large distance. Such profiles can change the distance to the last-scattering surface 
significantly, thereby allowing us to fit the CMB with a higher Hq. 

Finally we have shown that the model is consistent with the CMB dipole, if the position 
of the observer is at most about 0(10) Mpc/h away from the centre. On the other hand this 
would account for the recently measured large scale "bulk fiow" which seems inconsistent 
with ACDM. 

We have therefore demonstrated that there exist Void models based on LTB metrics, 
which may be a viable alternative to ACDM (although still with a rather low value of 
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the local Hubble parameter), and which ought to be tested and constrained by future 
observations. 

Note Added 

After the submission of the first version of the this manuscript two papers, [69] and [70], 
appeared in close succession which address similar issues. While a detailed comparison of 
the similarities and differences between the approaches and results of these three papers 
requires more time, we make here a few remarks which are already apparent. While we 
agree with [70] that SN+CMB+i?o rule out EdS based void models with simple profiles, we 
find that Curved void models can be consistent, although the value of Hq is still somewhat 
low compared to observations. In contrast [70] concludes that void models can be ruled 
out based on these three measurements. In our opinion there are four possible reasons for 
this disagreement: (i) We have an enlarged parameter space where we include background 
curvature. [70] has only looked at EdS voids, (ii) We consider "deep" voids where the 
central matter abundance can be less than 0.1. We did not find any convincing reason to 
exclude deep voids, which have also been considered previously in literature [31]. (iii) We 
considered a range of radial profiles. While not all profiles can produce a sufficiently large 
local Hubble parameter, some apparently can. (iv) Finally, we implemented an integrated 
MCMC approach. That is, we fit all the relevant data sets simultaneously with MCMC 
simulations. In contrast in [70] only SN and CMB were fit using MCMC and other con- 
straints were later applied by means of importance sampling, or compared a posteriori with 
the best fit parameters obtained from SN+CMB. When fitting only SN+CMB, the MCMC 
chain has no interest in going into areas in parameter space that embody a higher value for 
i7o,in or give a better fit to the BAO. In our work we find models that fit the CMB worse 
than ACDM, but fit the BAO better, leading to a net Xcs that is competitive with that of 
ACDM. When only considering CMB, these parameters are rejected. In other words, the 
authors of [70] have taken a prior on parameter space that reflects the favoured parameters 
by SN+CMB, and then fit this biased parameter space to Hq, BAO and other observ- 
ables. This does not necessarily lead to the same result as taking a completely flat prior on 
parameters and then simultaneously fitting CMB+BAO+SN+HST. This is probably the 
reason for having different values for -ffo,in in EdS voids quoted in this work and in [70]. 
Finally, we agree with [70], that void models may be in conffict with LSS measurements, 
although an analysis which considers perturbations of LTB would be needed to make firm 
conclusions, as we have stressed. 

The authors of [69] argue that if the void extends to the surface of last scattering, 
radiation can no longer be treated as a test field, and its effect on the metric should be 
included in the calculations. We acknowledge that our analysis and other analyses to date 
may be oversimplified in that respect. 
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A. Analytical Approximations for Conformal times and Redshifts 

LTB metrics form a class of exact solutions of General Relativity, which contains only dust, 
i.e. zero pressure source. They contain three arbitrary functions of the radial coordinate 
r: however one of them describes an inhomogeneous big bang time and we are not going 
to use it here. Another function can be eliminated redefining a new radial coordinate. So, 
we are left with one physically meaningful arbitrary function, k{r). The metric functions 
are known exactly, but only implicitly, and therefore one cannot find closed form analytic 
expressions for photon trajectories, t{r),z{r), in the LTB metric. Approximate analytical 
expressions exist in two perturbative regimes, either when (1) the inhomogeneities extend 
only a small distance compared to the Hubble radius, or when (2) the LTB space-time is 
only a small deviation from the EdS space-time. (1) corresponds to a Newtonian expansion 
which in our model is given by L ^ ''hoD while (2) is just the traditional cosmological 
perturbative expansion in small potentials and corresponds to having a small curvature 
function, k{r). 

The perturbative expansions have previously been used to analytically fit the super- 
novae data for small voids [14]. For our purpose, these analytical expressions provide 
important checks on our numerical simulations. Also, they provide physical insight into 
the various corrections to the redshift and distances in the presence of inhomogeneities. 

In the following, we use the definitions f = Mr and L = ML. 

A.l Next to leading Order New^tonian Approximation: 0{L?') & all orders in k 

In [14] Newtonian expansions were obtained for voids with small curvatures everywhere. 
Here we generalize this to profiles with general k{r) and also compute the corrections up 
to 0{L?') which is essential to computing the monopole correction due the presence of the 
LTB void. 

The general expressions for R and R' can be written as a power series in k{r): 



i?=^r-7V(l + X;i?2n7'V"F) , 



and 



3 



R! = 2a 



M 



oo 

1 



l + X](n+l)i?2„7'V"(rn' 



(A.l) 



(A.2) 
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The evolution equation for r is then given by (3.1) and (A. 2) 



dT _ a[l + ZT i?2n7'"T2-(ffe-)'] 



dr 



Vl + 2kf 



-a 



l + 5^ii2n7'"^'"(^^")' 



+ 0(f2). (A.3) 



From (A.3) we find the expression 

oo „ 

r = TE + n=TE-aY^ ^2n7'" / dr T2"(fA;")' , 

1 -^ 

The above integral can be evaluated as follows: 



(A.4) 



In= f dr r2"(fA;")' = r^'^rk'' - 2n f dr r^'^-^fk'' 



dr 
dr ' 



^2n^^n ^ 2na I dr r^^'-^fk'' 



r'^n-j^^ + 2nar2"-i f dr fk"^ 



1 + X^i22rn7""^""(^fc'")' 

1 

oo 
l + ^i?2™7'"V2-(r^'")' 



+ 0{f^) 
+ 0{f^) 



r-^n-j^n ^ 2nar2"-i 



/ df fk'^ + Y^ R2nn^'^T^"' [ dr rk'\rk'^)' 
Putting everything together we have 



(A.5) 



n = -a J^ i?2n7'" T^^fA;" + Inar^ 



2n-l 



/C)0 ^ 

df f/fc" + J]] i?2m7^"'ro2'" / dr fk^'irk" 



(A.6) 



Using the definitions of /„'s we can condense the above equation to 

Ti = -a (rfir^j^k) + 207^0 f dr rkfi{T^j^k) [l + /(ro^'fc) + rk'-f^T^fi{T^j^k)])jA.7) 

We cannot simplify the expression any further, but we should be able to compute it because 
both the functions /, /i are known, and so is the profile. Also note that at the boundary 
where k = 0, the first term do not contribute. 

We are now ready to look at the corrections to the redshift. From (3.2) we have 



dz 2a 



1 + Z T 



l + ^(n+l)/22„,7'V"(rA;")' 



dr . 



so that 



ln(l + z) = 2a 



fdr-+ [dr y(n+l)i?2„7'V"-nrA:")' 
Jo T Jq ^ 



{A.i 



(A.9) 



We can evaluate these integrals using similar techniques as for r: 
1 



dr - = dr — (l + — ) 

T Jq Te \ TEj Jo TE JQ 'q 



dr +o.ldr'^^^+Oif^) 
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Or, 



Jo = ^ ln(l + ze) + ^ I df ff{Tl^\) . 
For the second integral we have 



(A.IO) 



Jo -^0 L 1 

Putting everything together we have 

In (1 + z) = ln(l + ze) + ^ dr rf{T^j^k) + 2a V (n + l)i?2n7'"T2"- VA:'^ 
+ 2aV2"-2 V(n + l)(2n - I)ii2n7'" / 

^ JO 



+0(f3 



dr r/c" 



1 + Y, R2ml^"'r^"'{k"' + mffe'^-^A:') 



We can again rewrite the above expression in a more concise form: 



f2a^ /" 

1 + z = [1 + ze) exp < —2" / ( 

I '^0 ^0 



2a 



df rfir^Yk) + —[fir'Yk] + r'Ykfiir'Yk)] 



+ 



2a^ 



df r[3T^YkMT^Yk) + 2(To^7^A;)^/2(ro^7'A:) - fir^l'k)] 



Tq Jo 

X [l + /(roVfc) + fA;VT2/i(roVA:)]} 
Note that this is a useful 



(A.ll) 
(A.12) 



A. 2 0{k) & all orders in r 

In this subsection we will look at the opposite regime, when L can be arbitrarily large but 
the curvature is always small. Thus, we are going to keep only the linear order terms in k. 
The metric functions are thus approximated as 

R = ll^rr^l + R2ul) = ^7'^t2(1 + R2-f\^k) , 

vr 



and 



R'=^l^Al + R2l^r\rky], 



R'='^^^^\l + 2R2j'THrky 
9 r '- 



The evolution equation for r is then given by (3.1) 

-aM[l + R2-f^T^irky - k{Mrf] + ©(A:^) , (A.13) 



dT _ ^-f^M[l + R27^T^{rky] 



The above equation leads us to the iterative expression 



T = tq — ar — a 



R2I I dr T {rky — dr kr 
Jo Jo 



TE + n 



(A.14) 
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where ri now denotes the 0{k) corrections to the FLRW trajectory coming from the 
integrals. Let us look at the first integral: 

pf pf I pf 

/ = / df T (fk)' = T fk — df 2Tfk—- ~ Trjfk + a df 2TErk 
Jo Jo dr Jq 



T^rk + a dr 2Tork — a / dr 2r k . 
'0 Jo 



If we define the following functions 



then we have 



kn{f) = / drr^k(r) , 
'o 



/ = T^rk + 2aToki — 2a k2 . 



(A.15) 
(A.16) 



and therefore the approximate result 



Ti = —a [i?27 [T^rk + 2aToki — 2a /C2) — ^2] = —a [^27 (''"e^^ + 2aToki) — 6A;2/5] . 

(A.17) 
We emphasize that the above result is exact in f. 
We can now look at the redshift equation 



dz 



[1 + z)R' 



dz 



dr y/l + 2E l + z 
Thus we have 



2a [t~^ + 2R2-f'^T{kry - A;fV~^] df. 



ln(l + z) = 2a 



dr T + 2i?27 / dr T{kr)' —I dr kr t 
Jo Jo 



2_-l 



(A.18) 



Since we are only interested in 0{k) corrections, in the second and third integral we can 
replace r — )■ r^;, so that 

kf"^' 



ln(l + z) = 2a 



df r 

- + 2i?27' / dfTEikr)' 

T Jo 



dr 
TE 



2a[Ji + J2 - J3] • (A.19) 



J^ 



The first integral simplifies as 



T Jo 'TE 



teJ 



dr 



r -I 

df- ^ 



TE Jo "' r| 2a 



dr dr 



— ln{l + ZE)-^^+, ,_ „ 

2a aTE Jo cxte dr la 

2a aTE 

J_ln(i + ^^)-^L_:| + j3. 

2a aTE 2 



1 / X Ti 

ln(l + ze) 



ln(l + ZE 
df 



dr 



T\ 



aTE Jo TE 



Te 

[R2l^Tl{fk)' - kf^] 



drR2'j TE^rky 



dr 

— I 

'TE 



-kr 



To summarize, we have 



Ti 



ln(l + z) = ln(l + ze) - 2— + aJ2 . 

TE 



(A.20) 
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Now, J2 can be simplified as follows: 

J2 = 2R2-f^ / df TE{kr)' = 2R2-f^[TEkr + a df kf] = 2R2-f'^[TEkf + afci(f)] . 
Jo Jo 



Putting everything together we have the final formula for the redshift: 



1 + 2; = (1 + ZE)exp 



2n 

TE 



+ 2aR2-f'^TErk{r) + 2a^i?27^A;i(f) 



A. 3 Comparison bet^veen analytical approximations and numerics 



(A.21) 
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Figure 5: A comparison of numerical results with analytical approximations. Right: A comparison 
with Eq. (A. 12), keeping zb constant and varying Sq. For small zb there is perfect agreement. Left: 
A comparison with Eq. (A.21), keeping Sq constant and varying zb- Again, this time for small Sq, 
there is perfect agreement. 

Since at all times the metric functions and their derivatives are defined analytically 
by Eqs. (2.4,2.5,2.6), for any pair {r,t} we can calculate these functions up to very high 
accuracy by numerically inverting Eq. (2.6). The only less controllable calculation is the 
integration of the geodesic equations, Eqs. (3.1,3.2), backwards in time starting from an 
observer at r = 0, t = to- 

We have checked that the numerical code confirms the analytical approximations in 
the above-mentioned regimes. In figure 5 we see AT/T = 1 — TcMB.in/^CMB.cff, where 
^CMB.in = 2.726 K and TcMB,cflf is the value used for the calculation of the CMB power 
spectrum. On the left, a comparison with Eq. (A. 12) is made. For three values of zb, we 
show the difference as a function of increasing 6q. For tiny zb the difference is invisible, 
and for each of the three cases the ratio between the analytical and numerical result is a 
constant with respect to varying 5q, indicating that the discrepancy is dependent only on 
Zb- This was to be expected, as Eq. (A. 12) is at all orders in k. 

On the right in Figure 5 we compare the numerical results with those from Eq. (A.21). 
In this case we calculated AT/T for three fixed values of Jo ) as a function of a varying zb- 
Again, for small values of the expansion parameter, in this case small kmax (hence small 
60), the numerics and the analytics perfectly agree, and the relative (dis) agreement does 
not change with changing zb, as Eq. (A.21) is at all orders in r/L. 
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B. Determining the CMB temperature Tcs for an uncompensated void 

The exact relation between the observed CMB temperature and the CMB temperature to 
feed in camb is defined in equation (4.1), which we repeat here shghtly rewritten, 



rdec = reff(l + zdf)=robs(l + ^, 



ohsJ 



(B.l) 



Let eflf denote all quantities in the metric of the effective FLRW-observer, and obs denote 
all quantities in the metric of the real observer at the centre of the void. Then, 



feff 



1 -L .ydcc 



--T, 



1 -L ■ydcc 
^ ^ ^Cff 

°'^l + .,df l + Z^ffl + ^b. 



(B.2) 



Here z* denotes redshift z evaluated on the photon geodesic at time t*. The case where 



the metric outside of the void is exactly FLRW, we had 
the void that only asymptotically goes to FLRW, we have 

i+zfri+<bs 

and 



l+^doc 1+^.^ 






1. For 



(B.3) 






(B.4) 



with limt,_5.t^^^ e = 0. In the LTB metric, the geodesic equation describing the redshift can 
be written as 

dz S{r, t) 



Sir,t) 



(l + z), 



(B.5) 



such that 






1 




1 + 


Z{eS,ohs) (t) 




1 + 


Z*u 
obs 




1 + 


obs 



dt' 



'(eff,obs)l''(eflf,obs) 



it'),t') 



= (1 + 4 



to 1 + 2(eff,obs)(*') 'S'(eff,obs)(''(eff,obs)(i')) ^0 ' 

1 /■*''- dt' SoUrohs{t'),t') 



(B.6) 



obs J 



1 + <bs Jt, 1 + -2obs(i') SoUroh.{t'),t') 



dt' 



l + z. 



eff 



+ 



1 + 



l + <bs SoUrohs{t'),t') 

l + ZoUt')Sohs{roh.it'),t') 1 + Zes{t') Ses{t') ' I + Zes{t') S^sit' 



Sesit') 



'eff 



Sesit') 



l + z. 



eff 



1 -I- ^dec 
^ "•" ^eff 



tde 



dt' 



l + z*b3 >5obs(robs(i'),*') l + <ff 5cff(t') 



1 + ^obs(i') 5'obs(robs(i'), *') 1 + -^cff (i') 5'cff(t') 



(B.7) 
(B.8) 
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Therefore 



''dec 



dt' 



l + <bs Sobs{rohs{t'),t') l + <ff 5eff(t') 



1 + Zobs(t') Sohs{rohs{t'),t') 1 + Zeff(t') 5eff (f) 



(B.9) 
(B.IO) 



which obviously goes to zero for t* 
FLRW at r 



tdec- Since the void by construction converges to 



oo, tne Junction i+^^^^^(t,) s,Urobs{t'),t') 



1 + Z,s{t') Seff{t') 

cahy decreasing with time, and the integral can be estimated as 



must be monotoni- 



lel < 



''dec 



dt' 



u 



Zjf, X 



tjc X 



l + <bs S,Urohs{t'),t') 



l + <ff s,^{t'] 

1 + Zobs(t') 5obs(robs(t'), *') 1 + ^eff (to 5eff(iO 



1 + Z 



l + <bs ^obs(robs(tO,tO 

1 + Zobs(t') 5obs(robs(iO, t') 1 + Zeff(t') 5eff (f) 



eflf 



5'cff(i'] 



i'=t. 






(B.ll) 



We dropped r in ^^^tjt, as this ratio is independent of r. With equation (B.ll) we can 
integrate from the observer backwards in time until a time t* for which e is a few orders 



'-, such that up to high accuracy Teg 






It 



for voids embedded inside FLRW, at a time t^ for which the 



of magnitude smaller than 1 
is trivial to see that e 
photon is outside the void. 

It is also at this coordinate {r*^^, t*]-,s} that we can pretend the effective FLRW to start, 
in order to determine the cosmological parameters Q^ describing the FLRW cosmology of 
the CMB observer, as explained in section 4. 

C. The goodness of fit of a void model to a supernova observation 

C.l Likelihood of a multivalued function 

The quantities that are measured from each supernova in a supernova survey, are its redshift 
z and its angular diameter distance dx, or equivalently the luminosity distance di, = 
(1 + z)'^d\. The theoretical function that is fit to the data is hence dA{z). In an FLRW- 
universe, both redshift and time vary monotonically along a photon geodesic from the 
supernova to the observer, so there is no ambiguity when exchanging t, r and z as time 
parameters. For each value of z, there is a unique combination of {t, r} which determines 
the angular diameter distance to the observer as d^ = R{r, t). 

Inside the mass-compensated void however, the photon usually experiences a blueshift 
when crossing the dense shell that surrounds the void, as shown in the left panel in Fig. 6. 
This implies that even though coordinate distance r and cosmic time t vary monotonically 
along the geodesic, z no longer does so. The same value z is met up to three times along 
the same geodesic. As a consequence, for one value of redshift z, the theory predicts three 
possible outcomes for dA(z), as shown in the right panel of Fig. 6. 
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Q 



r (arbitrary units) 




1.1 1.2 1.3 



Figure 6: Le/t; Rcdshift z as a function of coordinate distance r along a photon geodesic pointing 
from an observer (r = 0) to a supernova, for a void with Sq — —0.8 and Zb=l. The photon 
experiences a blueshift when crossing the shell around z — 1. The definition of the coordinate size 
L of the void is such that the highest value of r for which z{r) — z^, is equal to r — L, hence the 
rightmost r in this figure for which z crosses the value of Zfc. Right: The angular diameter distance 
dAiz) for the same model, exemplifying the multiple theoretical predictions for the value of cIa 
corresponding to one value of z. 



Let us for consider a single data point for the moment. To construct the likelihood 
of an observation given a model M, we need to calculate the probability of measuring a 
value di when M is the underlying model, p{di\M). For a gaussian measurement error at 
measurement i, the probability of measuring a certain value di given a certain theoretical 
prediction Xi 



p{di\xi 



1 


r i{xi-d^y] 

[ 2 af \ 


V27rcr^ 



(C.l) 



Taking into account that the theory may predict a certain spread in Xj, such that for each 
Xi there is a probabiliy p{xi\M), the probability p{di\M) becomes 



p{di\M) = / dx pixi\M) 



1 



27ro"< 



exp 



1 {xi - dif 

2 a? 



(C.2) 



If the theory predicts only one possible outcome Xi, we simply have p{xi\M) = 5{xi — 
a^predicted)) such that Eq. (C.l) and Eq. (C.2) are identical, exchanging Xi and ^predicted- In 
our scenario however, for certain data points the theory predicts N values, hence 



N 



p{xi\M) = '^Cij6{xi - Xij) , 



(C.3) 



^Af 



normalized such that ^j=i Cij = 1- The indices ij denote the jth prediction at position 



i. For example, in our scenario 
the theory for a redshift z,, 



is the jth angular diameter distance predicted by 



J, and Cij is its corresponding weight. The likelihood of the 
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observation di then becomes 

p{d,\M) = Y, 



N 



Ci 



V 



i=i 



2Tra., 



exp 



2 



di)' 



(TT 



(C.4) 



For multiple data points, the likelihood £ of the data give the model is given by 

C = l[p{d,\M). (C.5) 

i 

For single valued predictions, carrying a subscript s for clarity, this reduces to the well 
known x^ statistic, 



Cs 



1 



(2vr)fnfa, 



exp 



1 ^ 

E 



{xi - di)' 



err 



Cie-'2^\ 



where Ci 



C-- 



In our scenario, the likelihood becomes. 



(27r)T 



N 

n 




1 (a;,; 



= Cie aXoff ^ 
where the last line defines the effective Xefr 



(C.6) 

(C.7) 
(C.8) 



-21n£ + 21nCi. 



C.2 Marginalization over Hq 

In the results of supernova observations, the normalization is unknown [25] . The normal- 
ization is set by choosing a value for Hq, such that a dimensionless Hod a is obtained. In a 
Markov Chain Monte Carlo analysis (MCMC), the normalization of the supernovae would 
influence the location of the best fit point in parameter space, thereby falsely interfering 
with the preferred values of Hq of other datasets in the same analysis, and falsely constrain- 
ing Hq around a certain value. Therefore we cannot ignore this effect during the analysis 
and marginalize over Hq after the MCMC has finished, but we have to marginalize over Hq 
already for the evaluation of the likelihood of each individual choice of parameters. This 
is well known [71, 35, 72], and is implemented in COSMOMC [35] by means of analytical 
marginalization. As in our scenario the likelihood function is different from the simple x^ 
function, also the marginalization changes. 

The analytical marginalization for the single valued case, with a Jeffreys prior on Hq, 
for gaussian errors in the distance modulus /i = 51og]^o [(1 + ^)^-f^o'^A(-z)] + 25 leads to 



£. 



dx Cc, 



:exp 



2^ ^ +2 U- ' 



-1 



\^ 



O". 



0"»- 



+ ln27r 



ilnVi 






(Tt 



(C.9) 
(CIO) 

(C.U) 
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with X = 5 logio -^0 • 

In the multivalued scenario, the product of the sum in the likelihood in Eq. (C.7) can be 
expanded, such that each term is exactly of the form of Eq. (C.6) and the marginalization 
integral can be performed for each term individually. However, if the theory predicts 
for example A^ values for M data points, this implies evaluating the relevant terms in 
Eq (C.ll) N'^ times, which becomes already a huge number for few datapoints. In practice 
this is impossible to do, and we resorted to a numerical evaluation of the marginalization 
integral. Taking into account that for the supernovae alone the likelihood is of the order 
— ln£ ~ 120, a rescaling had to be performed in the numerical integration. As it turns out, 
the numerical marginalization of the likelihood takes only O(10~^) s, where evaluating the 
N terms could take forever. 
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D. Figures and Tables 
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Figure 7: The marginalized ID posterior probabilities of the free parameters, on all of which 
we had a flat prior, when fitting to CMB-hBAO+SN+HST. ACDM in dashed, red, EdS Void in 
dashed-dotted, blue, and Curved Void in solid, black. Only fij. out is a derived parameter in the 
Void scenario, but since it is a free parameter for ACDM, we do plot it here for comparison. Note 
that our choice fcpjvot — 0.05 Mpc//i affects the central values of the scalar tilt and its running, 
which has no consequences. 
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Figure 8: The marginalized ID posterior probabilities of derived parameters, when fitting to 
CMB-hBAO-hSN-hHST. ACDM in dashed, red, EdS Void in dashed-dotted, blue, and Curved Void 
in solid, black. In ACDM we have ifo.in = ^o out- The variable ifo.in is fit to the HST value. 
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Table 11: The best fit parameters (each left column) and marginalized posterior probabilities 
(each right colmTin) with 95% C.L. errors of all free parameters (top ten parameters) and derived 
parameters (from the eleventh parameter down) , for the three models at stake when simultaneously 
fitting CMB + BAO + SN + HST. In the Curved Void-model, the parameters ftk is actually a 
derived parameter, as we take a flat prior on ^dm.out in that case (see Table 1). Let us again 
emphasize that our pivot scale fcpjvot = 0-05 Mpc / /i, which leads to the lower value of rig, since we 
allow for a running. 



59 




0.02 0.025 

a, h- 



CMB + SN + LSS + HST 




l\ 



0.1 0.15 

a, h^ 

dm. uul 



0.2 



II /^ 



40 



60 




0.04).0ffl.080.10.123.14 

t 





-0.5 -0.4 -0.3 




0.9 1 1.1 1.2 1.3 





3.2 3.4 

log[10'"A] 




/\ 



-0.2 -0.1 



Figure 9: The marginalized ID posterior probabilities of the free parameters, on all of which 
we had a flat prior, when fitting to CMB+SN+SDSS-^HST. ACDM in dashed, red, EdS Void in 
dashcd-dotted, blue, and Cm^ved Void in solid, black. Only flk.out is a derived parameter in the 
Void scenario, but since it is a free parameter for ACDM, we do plot it here for comparison. Note 
that our choice fcpivot — 0.05 Mpc / h affects the central values of the scalar tilt and its running, 
which has no consequences. 
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Figure 10: The marginalized ID posterior probabilities of derived parameters, when fitting to 
CMB-hSN-l-SDSS+HST. ACDM in dashed, red, EdS Void in dashed-dotted, blue, and Curved Void 
in solid, black. In ACDM we have i/o.in = -f^o.out- The variable i^o.in is fit to the HST value. 



60 



CMB+SN+SDSS+HST 



Parameter 


ACDM 


Cur 


ved Void 


EdS Void 


^b,oath^ 


2.21 X 10-^ 


2.12+|;;;5 X 10-2 


2.24 X 10-2 


2.i9«:;» X 10-2 


1.90 X 10-2 


1.92+[!;;g X 10-2 


f'^dm.oath 


0.111 


U.i2/_m,i7 


0.148 


r, 1 ^r +0.020 
0-145_o,oi7 


0.186 


n is«+""ii 

0.i»S_(,|,„(j 


-ff(),out |kmB-' Mpc-'] 


56.7 


60.0+Ii 


38.3 


38-0+2:9 


45.2 


45.6«J 


r 


8.44 X 10-^ 


8-86+2:76 X 10"^ 


8.13 X 10-2 


8.40+2;!^:; X 10-2 


7.57 X 10-2 


8-43^2:89 X 10-2 


i^fc,out 


-3.50 X 10-2 


-l.45l;li^ X 10"' 


-0.162 


-0.156«;!!^? 


- 


- 


"s 


0.970 


1 n3+°-i3 


1.06 


i-07lg:;i! 


1.25 


1 nr+0.07 

1-^0-0.07 


"s 


-1.08 X 10-2 


-4.30«;:i^ X 10-' 


-6.09 X 10-2 


-6.36lf ?i X 10-2 


-0.158 


-n 1 c;fi+co3i 


log[10io^] 


3.20 


o 01 +0.08 
•'■^-'--0.08 


3.25 


o , j+0.08 


3.26 


n no+0.08 

-^■^8-0.08 


ZB 


- 


- 


2.90 


1 q-l+1.09 


1.16 


1 49+l'll 
-'-■^''-0.76 


So 


- 


- 


-0.465 


n 4fi7+''"6."i 
-0.4b('_(,(|gj 


-0.449 


n 4c:-i +0.068 


fU 


0.622 


n 61 4+""57 


- 




- 




' '7n,out 


0.413 


0.400+°:°?J 


1.16 


i-i6+g:g? 


1.00 


1.00 


-f^O.in [kms-' Mpc-'l 


56.7 


60.0+Ii 


45.6 


45.1+1^ 


52.7 


r,o n+2.2 
^'^■'^-2.0 


T2.72/Teff 


1.00 


1.00 


1.01 


0.996«:gJ^ 


0.990 


n qqo+0.008 


*o [Gyr] 


15.3 


14.6«:2 


16.5 


16.7+- 


14.4 


14.3;«i 


-ffo.eff [km ■"' Mpc-'] 


56.7 


60.0+^:? 


36.4 


37.1«i 


45.2 


45.2+S:? 


nt,esfi' 


2.21 X 10-2 


2.12l°;;5 X 10-2 


2.03 X 10-2 


2.09+1]:;^ X 10-2 


1.90 X 10-2 


1.89;U;™ X 10-2 


f2d,„,eff/l.' 


0.111 


0.122tS:S? 


0.134 


n 1 Qs+0.019 
0-J-3S-0.017 


0.186 


n 1 oR+0.005 
0-18o_o.oo5 


tJt.eff 


-3.50 X 10-2 


-1.45l|^^ X 10-2 


-0.168 


-0.i58t!!;!j?,^ 


- 


- 


n„,eff 


0.413 


U.400_Qf,yQ 


1.17 


i-i6^S:S? 


1.00 


1.00 


i\Jl^ 


2.21 X 10-2 


2.12+g:;^ X 10-2 


1.20 X 10-2 


i.i7i|!:;;i X 10-2 


1.05 X 10-2 


1.05;J!;;2 X 10-2 


i1dm,inh? 


0.111 


r, -|nn+0.019 

U.i22_o(,i7 


7.93 X 10-2 


7-70i;:25 X 10-2 


0.102 


0.103«;||;2 


f2;t,in 


-3.50 X 10-2 


-1.45l|^!! X 10-2 


0.560 


c;Rq+0.067 


0.594 


0.596tSi;|!« 


n„,in 


0.413 


n 400+'^°''i 

U.400_op7Q 


0.440 


r, ^07+0.073 
O-^' -0.067 


0.406 


0.404i;«^2 



Table 12: The best fit parameters (each left column) and marginalized posterior probabilities 
(each right column) with 95% C.L. errors of all free parameters (top ten parameters) and derived 
parameters (from the eleventh parameter down), for the three models at stake when simultaneously 
fitting CMB + SN + SDSS + HST. In the C-model, the parameters flj. is actually a derived 
parameter, as we take a flat prior on ri^m, out in that case. Let us again emphasize that our pivot 
scale fcpivot = 0.05 Mpc / h, which leads to the lower value of Ug, since we allow for a running. 
The most notable difference with Table 11 is the smaller value of the density contrast ^o- Since 
we use FLRW perturbation theory for the LSS inside the void, a fit to the matter power spectrum 
has considerable weight on the value of VL^i^ inside the void. Note that this approximation is not 
justified. 
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Table 13: The best fit parameters of all free parameters (top twelve parameters) and derived 
parameters (from the thirteenth parameter down), for the ACDM and all considered curvature 
profiles when simultaneously fitting CMB + BAO + SN + HST. The parameter 5o(fcmax-l-fcmax,2, ^o) 
is the actual density contrast between the centre of the void and the outer FLRW, where in the 
case of Profile E the outer FLRW is never reached. The parameter ^(fcmaxjio) gives the density 
contrast that would be present if fcmax,2 — 0, hence defining both /cmax and fcmax.2 for a given 
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Figure 11: Plots of all the datasets and the best fit models when fitting to CMB + BAO + SN + 
HST. The top panel shows SN data and predictions normalized to a reference model with 51„i — 0.3 
and fife =0.7 (OCDM). The second panel from the top shows CMB data and predicitions normlized 
to the ACDM best fit model for this combination of datasets. The bottom left panel shows the 
BAO, the bottom right panel i?o,in- In all plots ACDM is dashed green, the Curved Void-model 
dotted blue, and the EdS Void-model fine dotted magenta. 
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Figure 12: Plots of all the datasets and the best fit models when fitting to CMB + SDSS + SN + 
HST. The top panel shows SN data and predictions normalized to a reference model with i7„j ~ 0.3 
and fife =0.7 (OCDM). The second panel from the top shows CMB data and predicitions normlized 
to the ACDM best fit model for this combination of datasets. The bottom left panel shows -ffo.in, 
and the bottom right panel shows the galaxy power spectrum, as well as the contribution to the 
total x^ for each separate data point, in a bar chart. In all plots ACDM is dashed green, the Curved 
Void-model dotted blue, and the EdS Void-model fine dotted magenta. 
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